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Abstract. We present new or improved methods for calculating NLTE, line-blanketed model atmospheres for hot stars with 
winds (spectral types A to O), with particular emphasis on a fast perfor mance. These methods hav e been implemented into a 
previous, more simple version of the model atmosphere code Fastwlnd iSantolava-Rev et alll997^ and allow to spectroscop- 
ically analyze rather large samples of massive stars in a reasonable time-scale, using state-of-the-art physics. Although this 
updated version of the code has already been used in a number of recent investigations, the corresponding methods have not 
been explained in detail so far, and no rigorous comparison with results from alternative codes has been performed. This paper 
intends to address both topics. 

In particular, we describe our (partly approximate) approach to solve the equations of statistical equilibrium for those elements 
which are primarily responsible for line-blocking and blanketing, as well as an approximate treatment of the line-blocking 
itself, which is based on a simple statistical approach using suitable means for line opacities and emissivities. Both methods are 
validated by specific tests. Furthermore, we comment on our implementation of a consistent temperature structure. 
In the second part, we concentrate on a detailed comparison with results from those two codes which have b een us ed in 
alternative spectroscopical investigations, namely cmfgen ( Hillier & Miller 1998) and WM-Basic I Pauldrach et all EoOlh . All 
three codes predict almost identical temperature structures and fluxes for A > 400 A, whereas at lower wavelengths a number 
of discrepancies are found. Particularly in the Hen continua, where fluxes and corresponding numbers of ionizing photons react 
extremely sensitively to subtle differences in the models, we consider any uncritical use of these quantities (e.g., in the context 
of nebula diagnostics) as being dangerous. Optical H/He lines as synthesized by fastwind are compared with results from 
CMFGEN, obtaining a remarkable coincidence, except for the Hei singlets in the temperature range between 36,000 to 41,000 K 
for dwarfs and between 31,000 to 35,000 K for supergiants, where cmfgen predicts much weaker lines. Consequences due to 
these discrepancies are discussed. 

Finally, suggestions are presented how to adequately parameterize model-grids for hot stars with winds, with only one additional 
parameter compared to standard grids from plane-parallel, hydrostatic models. 

Key words. Methods: numerical ~ Line: formation - Stars: atmospheres - Stars: early-type ~ Stars: mass-loss 



1. Introduction 

During the last years, the quantitative spectroscopy of mas- 
sive stars with winds has made enormous progress due to 
the development of NLTE (non-local thermodynamic equi- 
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librium) atmosphere codes which allow for the treatment of 
metal-line blocking and blanketing. With respect to both spec- 
tral range (from the extreme ultraviolet, EUV, to the in- 
frared, IR) and metallicity of the analyzed objects (from SMC- 
abundances to Galactic center stars), a wide range in pa- 
rameters can now be covered. Presently, five different codes 
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are in use which have been developed for specific objec- 
tives, but due to constant improvements they can be applied 
in ot her contexts as well. In particular, these codes are cm- 
fgen (lHimer&Mil lei<'l998l). the "Potsdam-group" code de- 
veloped b^_W;R;_^amann and collaborators (for a status re- 
port, see lGrafener et al. 2002), th e "multi-p urpose model atmo- 
sphere code" ph oenix (Hausch ildt & Baroa.1999.) . WM-Basic 
jPauldrach et al. |^01) and fas twdsto , which wil l be described 
here (see also Santolava-Rev et alJ 1 19971 and iHerrero et alJ 
l2002l for previous versions). 

The first three of these codes are the most "exact" ones, 
since all lines (including those from iron-group elements) are 
treated in the comoving frame (CMF), which of course is a very 
time-consuming task. Moreover, since the first two of these 
codes have originally been designed for the analysis of the 
very dense winds from Wolf-Rayet stars, the treatment of the 
photospheric density stratification is approximative (constant 
photospheric scale-height). For several analyses this problem 
has been cured by "coupling" cmfgen with the plane-parallel, 
hydr ostatic code tlusty developed by Hubenv & Lanz ( 1995) 
(e.g.. lBouret et alJ2003l) . 

The multi-purpose code phoenix is mainly used for the anal- 
ysis of supernovae and (very) cool dwarfs, but also a small 
number of hotter objec ts have been considered , e.g., the A- 
type supergiant Deneb jAufdenberg et alJl2002 l). Due to this 
small number a detailed comparison with corresponding results 
is presently not possible, and, therefore, we will defer this im- 
portant task until more material becomes available. 

In contrast to all other codes which use a pre-described 
mass-loss rate and velocity field for the wind structure, the 
model atmospheres from WM-Basic are calculated by actually 
solving the hydrodynamical equations (with the radiative line- 
pressure being_a22roximate^ con- 
cept, cf. lCastor et al.ll975tlPauldrach et al.ll986 ) deep into the 
photosphere. Thus, this code provides a more realistic stratifi- 
cation of density and velocity, particularly in the transonic re- 
gion (with the disadvantage that the slope of the velocity field 
cannot be manipulated if the wind does not behave as theoret- 
ically predicted). Since WM-Basic aims mainly at the predic- 
tion of EUV/UV fluxes and profiles, the bound-bound radiative 
rates are calculated using the Sobolev approximation (includ- 
ing continuum interactions), which yields "almost" exact re- 
sults except for those lines which are formed in the transonic re- 
gion (e.g., Santolava-Rev et al. 1997). Moreover, line-blocking 
is treated in an effective way (by means of opacity sampling 
throughout a first iteration cycle, and "exactly" in the final iter- 
ations), so that the computational time is significantly reduced 
compared to the former three codes. 

FASTWiND, finally, has been designed to cope with optical 
and IR spectroscopy of "normal" stars with Tg^ ^ 8,500 K', 
i.e., OBA-stars of all luminosity classes and wind strengths. 

Since the parameter space investigated for the analysis 
of one object alone is large, comprising the simultaneous 
derivation of effective temperature Tetf , gravity log^, wind- 
strength parameter Q = M/(/?*Vco)'^ (cf. Sect. |9}, veloc- 



ity field parameter p, individual abundances (most important: 
helium-abundance Yue ) and also global background metallic- 
ity z, much computational effort is needed to calculate the 
vast amount of necessary models. This is one of the reasons 
why the samples which have been analyzed so far by both 
CMFGEN and WM-Basic are not particularly large^ , comprisin 



typically five to sev en objects per analysis (e.g.. Sillier et alJ 
|2003; Bour etet all2 003: Martins et al. 2004 for recent cmfgen- 
analys es and Fullerton et al.l I2000t iBianchi & Garcial l2002t 
lOarcia & BianchL2004 for recent WM-Basic analyses). 

The reader may note that although the number of fit- 
parameters gets smaller when the wind-strength becomes neg- 
ligible, a certain difference between the results from "wind- 
codes" and plane-parallel, hydrostatic model atmospheres still 
remains: More or less independent of the actual mass-loss rate, 
there will always be an enhanced probability of photon es- 
cape from lines in regions close to the sonic point and above, 
if a super-sonic velocity field is present. A prime example il- 
luminating the consequences of this enhanced escape is th e 
Hen ground-state depopulation in O-stars (Gableret al] ll989l) . 
even though it is diminished by line-blocking effects compared 
to the original case studied by means of pure H/He atmospheres 
(see also Sect. 14. 7> . 

With the advent of new telescopes and multi-object spec- 
trographs, the number of objects which can be observed during 
one run has significantly increased (e.g., flames attached to the 
VLT allows for an observation of roughly 120 objects in paral- 
lel). An analysis of those samples will definitely result in more 
reliable parameters due to more extensive statistics but remains 
prohibitive unless the available codes are considerably fast. 

This was and still is the motivation which has driven the de- 
velopment of FASTWIND. We have always considered a fast per- 
formance to be of highest priority. The required computational 
efficiency is obtained by applying appropriate physical approx- 
imations to processes where high accuracy is not needed (re- 
garding the objective of the analysis - optical/IR lines), in par- 
ticular concerning the treatment of the metal-line background 
opacities. 

Meanwhile, a number of analyses have been performed 
with our present version of fastwind, with significant sam- 
ple sizes, of the ord e r of 10 to 40 stars per sample (e.g., 
Urbaneia et ail l2003t iTrundle et alJ l2004t lUrbaneid l2004 
alJ l2004t iMassev et all |2004 l2005 l). Although 
the code has been carefully tested and first comparisons 
with results from c mfgen and tlusty have been published 
(lHerreroetalJ l2002'). a detailed description of the code and 
an extensive comparison have not been presented so far 
Particularly the latter task is extremely important, because oth- 
erwise it is almost impossible to compare the results from anal- 
yses performed using different codes and to draw appropriate 
conclusions. An example of this difficulty is the discrepancy 
in stellar parameters if results from optical and UV analyses 
are compared. Typically, UV-spectroscopy seems to result in 
lower values for Teff than a corresponding optical analysis, e.g.. 



' i.e, molecules do not play any role and hydrogen remains fairly 
ionized. 



- From here on, we will concentrate on the latter two codes because 
of our objective of analyzing "normal" stars, whereas the "Potsdam"- 
code has mainly been used to analyze WR-stars. 
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iMassev et al](l2005l) . Unless the different codes have been care- 
fully compared, no one can be sure whether this is a problem 
related to either inadequate physics or certain inconsistencies 
within the codes. 

This paper intends to answer part of these questions and is 
organized as follows: In Sect. |2 we give a quick overview of 
the basic philosophy of the code, and in Sect.|3lwe describe the 
atomic data used as well as our treatment of metallicity regard- 
ing the flux-blocking background elements. Sects. I^andlsjgive 
a detailed description of our approach to obtain the fast perfor- 
mance desired: Sect.|3]details on the approximate NLTE solu- 
tion for the background elements (which is applied if no con- 
sistent temperature structure is aimed at), and Sect.|5]describes 
our present method to tackle the problem of line-blocking. Both 
sections include important tests which have convinced our- 
selves of the validity of our approach, particularly after a com- 
parison with results from WM-Basic. Sect.|6lcovers the problem 
of level inversions and how to deal with them, and Sect.0com- 
prises the calculation of a consistent temperature structure. In 
Sect. |8] a detailed comparison with results from a grid of cmf- 
GEN models'* is performed, and Sect.|9lsuggests how to param- 
eterize model-grids adequately and reports on first progress. In 
Sect. [TO] finally, we present our summary and an outlook re- 
garding future work. 

2. Basic philosophy of the code 

In the following, we will summarize the basic features of 
FASTWiND, before we describe in detail the methods used. The 
first version of the code (unblocked atmosphe r e/line formation) 
has been introduced bv ISantolava-Rev et al.l ( 1 19971 hereafter 
Paper I), and has been significantly improved meanwhile. Let 
us first mention that we distinguish between two groups of el- 
ements, namely the so-called explicit ones and the background 
elements. 

The explicit elements (mainly H, He, but also C, N, O, Si, 
Mg in the B-star range, see below) are those which are used 
as diagnostic tools and are treated with high precision, i.e., by 
detailed atomic models and by means of CMF transport for the 
bound-bound transitions. In order to allow for a high degree 
of flexibility and to make use of any improvements in atomic 
physics calculations, the code is atomic data driven with respect 
to these ions, as explained in Paper I: the atomic models, all 
necessary data and the information on how to use these data are 
containe d in a user supplied file (i n the so-called detail input 
form, cf. iButler & Giddingslll985h whereas the code itself is 
independent of any specific data. 

The background ions, on the other hand, are those allowing 
for the effects of line-blocking/blanketing. The corresponding 
data originate from Pauldrach et al. ( 1998, 20()1J and are used 
as provided, i.e., in a certain, fixed form. 

FASTWIND follows the concept of "unified model atmo- 
spheres" (i.e., a smooth transition from a pseudo-hydrostatic 
photosphere to the wind) along with an appropriate treatment 
of line-broadening (Stark, pressure-) which is a prerequisite 
for the analysis of O-stars of different luminosity classes cov- 



ering a variety of wind densities. Particularly and as already 
described in Paper I, the photospheric density consistently ac- 
counts for the temperature stratification and the actual radia- 
tion pressure, now by including both the explicit and the back- 
ground elements. 

The corresponding occupation numbers and opacities (of 
the background-elements) can be derived in two alternative 
ways: 

a) in those cases, where the temperature stratification is cal- 
culated by means of NLTE Hopf parameters (see be- 
low), we apply an approximate NLTE solution for all 
background elements foll owin g the p ri ncipal ph il osoph y 
developed by [A bbott & I^ ( 1 1985b. ISchmutj ( Il99ll) . 
'Schae rer & Schmutzi (i 19941) and lPuls et alJ ilOOd) . where 
important features have now been improved (cf. Sect. l4}. 
Particularly, the equations of approximate ionization equi- 
librium have been re-formulated to account for the actual 
radiation field as a function of depth and frequency, and a 
consistent iteration scheme regarding the coupling of the 
rate equations and the radiation field has been established 
to avoid the well-known convergence problems of a pure 
Lambda Iteration ("Sect. l4~6l . 

b) In the other case, when the T-stratification shall be calcu- 
lated from first principles, the complete set of rate equations 
is solved almost "exactly" for the most abundant back- 
ground elements (C, N, O, Ne, Mg, Si, S, Ar, Fe, Ni, if not 
included as explicit ions), employing the Sobolev approx- 
imation for the net radiative rates (with actual illuminat- 
ing radiation field). The remaining background elements, 
on the other hand, remain to be treated by the approxima- 
tion as outhned in a). 

In order to account for the effects of line-blocking, we use suit- 
able means for the line opacities, averaged over a frequency 
interval of the order of 1,000. . . 1,500 km s and appropriate 
emissivities (Sect.lS}- 

Finally, the temperature stratification can be calculated in 
two different ways. As long as one is exclusively interested 
in an optical analysis, the concept of NLTE-Hopf parameters 
(cf. Paper I) is still sufficient, if the background elements are 
accounted for in a consistent way, i.e., have been included in 
the particular models from which these parameters are derived. 
Since this method is flux-conservative, the correct amount of 
line-blanketing is "automatically" obtained. Note that for op- 
tical depths tr ^ 0.01 a lower cut-off temperature is defined, 
typically at T^in = 0.6reff. 

Alternatively, the new version of fastwind allows for 
the calculation of a consistent"^ temperature, utilizing a flux- 
correction method in the lower atmosphere and the thermal 
balance of electro ns in the outer one (Sect. 0. As has been 
discussed, e.g., bv lKubat et alJ ( ll999l) . the latter method is ad- 
vantageous compared to exploiting the condition of radiative 
equilibrium in those regions where the radiation field becomes 
almost independent on Tg. Particularly for the IR-spectroscopy, 
such a consistent T-stratification is important, since the IR is 



as recently calculated bv lLenorzer et alj i2004) 



* Note, however, that non-radiative heating processes might be of 
importance, e.g., due to shocks. 
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formed above the stellar photosphere in most cases and de- 
pends critically on the run of in those regions, where our 
first method is no longer applicable. 



3. Atomic Data and Metallicity 

Explicit elements. In order to obtain reliable results also in 
the IR, we have significantly updated our H- and He-models 
compared to those described in Paper I. Our present H and 
Hen models consist of 20 levels each (vs. 10 and 14 in the 
previous version, respectively), and Hei includes levels until 
n = 10, where levels with n - 8. . .10 have been packed (previ- 
ous version: 8 levels, packed from 5. . . 8). Fur ther informatio n 
concerning cross-sections etc. can be found in Jokuthv '("2002). 
Present atomic models for metals have been accumulated from 
different sources, mainly with respect to an analysis of B-stars, 
i.e., for ionization stages n and m, except for Mg (i,n) and Si 
(ii. III, iv). Informatio n on our Si atomic model can be found in 
iTrundle et all ( l2004 . and on the oth er metals incorporated so 
far (C, N, O, Mg) in lUrbaneial l Eool . 

Background elemen ts. The atomic data for back ground ele- 
ments originate from lPauldrach etal] lll998lEoOll) . who have 
given a detailed description on the various approaches and 
sources. These data comprise the elements from hydrogen to 
zinc (except Li, Be, B and Sc which are too rare to affect the 
background opacity) with ionization stages up to viii. The num- 
ber of connecting lines (lower and upper level present in the 
rate equations) is of the order of 30,000, and the number of 
lines where only the lower level is present is 4.2 ■ 10^. The for- 
mer group of lines is used to solve the rate equations, whereas 
the latter is used to derive the metal-hne background opacities 
(cf. Sect In addition to bound-free cross-sections and gf- 
values, there is also detailed information about the collision- 
strengths for the most important coUisional bound-bound tran- 
sitions in each ion. 

Metallicity. The abundances of the backg round elements are 
taken from the solar values provided by iGrevesse & Sauvall 
lll998i and references therein)^. For different "global" metal- 
licities, z - ZIZq, these abundances are scaled proportionally 
with respect to mass ratios, e.g., by 0.2 for the SMC and by 0.5 
for the LMC (although these values are certainly disputable, 
e.g. Ma ssev et al. 2004 and references therein). 

A particular problem (independent on the actual value of z) 
appears in those cases when the He/H ratio becomes non-solar. 
In this case, we retain the specific relative mass fractions of 
the other elements, which of course has a significant effect on 
the number ratios. Although this procedure is not quite right, it 
preserves at least the overall mass fraction of the metals, par- 
ticularly the unprocessed iron group elements, which are most 
important for the line-blocking. Further comments on the valid- 
ity of this procedure have been given bv .Massev et al. ( .2004j) . 
Let us briefly mention a comparison to evolutionary calcula- 

' Of course, the user is free to change these numbers. 



tions from lSchaerer et al I d 19931) performed by P. Massey (priv. 
comm.). 

For the 120 Mq track at Z = 0.008 (roughly the LMC metal- 
licity), Z stays essentially unchanged in the core until the end 
of core H burning, even though the mass fraction of C and N 
are going up while O is going down: at a number ratio Fee - 
2 (i.e., the mass ratio Y has changed from 0.265 to 0.892), the 
value for Z has changed insignificantly from 0.0080 to 0.0077, 
and even more interestingly, the mass fraction of the sum of C, 
N, O, and Ne has essentially changed in the same way (0.0075 
to 0.0070), even though the actual mass fraction of N has more 
than doubled. 

4. Background elements: approximate NLTE 
occupation numbers 

In order to save significant computational effbrt, the occupa- 
tion numbers of the background elements are calculated by 
means of an approximate solution of the NLTE rate equa- 
tions. Such an approach has been successfully applied in a va- 
riety of stellar atmosphere calculations, e.g., to derive the ra- 
diative acceleration of hot star winds dAbbott & LucvllT985[ 
ILucv & Abbott I993I. and for the spe ctroscopy of hot stars 
(iSchmutzl | 199U iSc haerer & Schmutz 1994) and Supernova 
remnants (*Mazzali & Lucv..l993; .Lucv..l999.: MazzalL.2O()0l) . 
Puis et al. (2000) have used this method for an examination 
of the line-statistics in hot star winds, by closely following a 
procedure discussed by lSpringmannl(ll997l) which in turn goes 
back to unpublished notes by L. Lucy. 

One might argue that such an approximate approach can 
barely handle all the complications arising from sophisticated 
NLTE effects. However, in the following we will show that the 
approximate treatment is able to match "exact" NLTE calcula- 
tions to an astonishingly high degree, at least if some modifica- 
tions are applied to the original approach. Moreover, the calcu- 
lated occupation numbers will not be used to synthesize line- 
spectra, but serve "only" as lower levels for the line-opacities 
involved in the blocking calculations. 

Actually, the major weakness of the original approach is 
the assumption of a radiation field with frequency independent 
radiation temperatures T^^i- Since solely the difference in ra- 
diation temperatures at strong ionization edges is responsible 
for a number of important effects, we have improved upon this 
simplifiction by using consistent radiation temperatures (taken 
from the solution of the equations of radiative transfer). As we 
will see in the following, this principally minor modification 
requires a number of additional considerations. 

4.1. Selection of levels 

One of the major ingredients entering the approximate solu- 
tion of the rate equations is a careful selection of participat- 
ing atomic levels. In agreement with the argumentation by 
lAbbott & Lucvl ll 19851) only the following levels are used: 

- the ground-state level 

- all meta-stable levels (from equal and different spin sys- 
tems), denoted by "M" 
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- all excited levels which are coupled to the ground-state via 
one single permitted transition where this transition is the 
strongest among all possible downward transitions; in the 
following denoted as subordinate levels "N" 

- all excited levels coupled to one of the meta-stable levels 
OT e M in a similar way (subordinate levels "S")- 

In the above definition, the term "strongest" refers to the 
Einstein-coefficients Ay,. All other levels are neglected, since 
their population is usually too low to be of importance and can- 
not be approximated by simple methods. 



4.2. Ionization equilibrium 

In order to allow for a fast and clearly structured algorithm, 
we allow only for ionizations to and recombinations from 
the ground-state of the next higher ion, even if this is not 
the case in reality. Due to this restriction and by summing 
over all line-processes an "exact" rate equation connecting 
two neighboring ions is derived which exclusively consists of 
ionization/recombination processes. In the following, we will 
further neglect any collisional ionization/recombination pro- 
cesses, which is legitimate in the context considered here, 
namely in the NLTE-controlled atmospheric regime of hot 
stars. (In the lowermost, LTE dominated part of the atmo- 
sphere, tr > 2/3, we approximate the occupation numbers a 
priori by LTE conditions). 

At first, let us consider an ion with only one spin system, 
e.g., a hydrogenic one. In this case, the ionization equilibrium 
becomes 



(1) 



with n, the occupation numbers of the lower ionization stage, 
riic the (ground-state!) occupation number of the higher ion, the 
asterisks denoting LTE- conditions (at a ctual electron density, 
(n,/nj* = tiMJ), cf. lMihalasl[l97l Sect. 5) and ioniza- 
tion/recombination rate coefficients 



Rn 



AnJy 
hv 



-fl(v)dv 



a{v) dv. 



(2) 



(3) 



Jy is the mean intensity, a{y) the ionization cross-section and all 
other symbols have their usual meaning. Once more, within our 
above approximation (ionization to ground state only), Eq.^is 
"exact" and does not depend on any assumption concerning the 
bound-bound processes (radiative or collisional; optically thick 
or thin) since the corresponding rates drop out after summation. 

By introducing the recombination coefficient a, defined in 
the conventional way. 



n^neUi = n.R^i, 

the ionization equilibrium can be reformulated 

^ = n^ne ^ Ui, 



(4) 



(5) 



and we extract all quantities referring to the ground-state of the 
lower ion. 



n\R\K\l + 

n\RiK 



(6) 



Finally, inserting the ground-state recombination coefficient a\ 
(cf. Eq. 0} on the rhs, we obtain the ionization equilibrium ex- 
pressed as the ratio of two neighboring ground-states. 











1 RkI 





1 + 



(7) 



Note that this ratio depends on the radiation field, the actual 
electron density and temperature and on the excitation within 
the lower ion, which will be discussed in the next subsection. 
Note also that all what follows is "only" a simplification of this 
equation. 

So far, our derivation and the above result are identical 
to previous versions of the approximate approach. From now 
on, however, we will include the frequency dependence of the 
radiation field. To this end, we descri be the ionizat ion cross- 
sections by the Seaton-approximation ( ISeatODlll958l) . which is 
not too bad for most ions. 



a(v) 



(8) 



Writing the mean intensity as Jv(r) - W(r)Bv(Tad(v, r)) with 
dilution factor W and neglecting the stimulated emission in 
the recombination integral (valid for all important ionization 
edges), we obtain (radial dependence of all quantities sup- 
pressed in the following) 



Ri, 



IT 



aiT{x,x,l3, s) 



with Xr i = hvilkT^^i, jCe,i - hvi/kTf. and 

r(x;/3, s)^/3x'T(3-s,x) + (l -/?) ji:'+T(2 - x). 



(9) 

(10) 

(11) 



We have assumed Tr^i{v) -: T^ i to be constant over the deci- 
sive range of the ionizing continuum v ^ v,, since only those 
frequencies close to the edge are relevant. In other words, each 
transition is described by a unique radiation temperature. In the 
above equation, the incomplete Gamma-function T{a,x) has 
been generalized to include also negative parameters, a < 0. 
The ratio of ground-state ionization/recombination rate coeffi- 
cients is thus given by 



^.1 I Te / r(x^x,i3,sy 



(12) 



i.e., is independent of the actual value of the cross-section at 
the threshold, a,-. Although this expression is rather simple, it 
requires the somewhat time-consuming evaluation of the in- 
complete Gamma-functions. To keep things as fast as possible, 
we generally use the parameter set 06 = 1 , i = 2) instead of the 
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Fig. 1. Ratio of ionization to recombination rate coefficients: relative error between "exact" ratios rEq. ll2> and approximate ones 
(Eq. \T5\ with /? = 1 and s = 2) as a function of T^^d/Te, for different combinations of (yS, s ). The error decreases for even higher 
ionization energies. 



actual parameters which results in a particularly simple func- 
tion T, 

'Fix;l,2)^x^exp(-x). (13) 

Note that these parameters do not correspond to the hydrogenic 
cross-section, which would be described hy s - 3. Using this 
parameter set, the ionization/recombination rates simplify to 



hvi 













(14) 



(15) 



SnIkTA J 
and the ratio R\kIRk\ becomes 

Rik 1 

— - W — '— exp 

We have convinced ourselves that this approximation leads to 
acceptable errors of the order of 10%, cf. Fig.^ Furthermore, 
we define the following quantities, where ^ is just the ratio of 
ground-state to total recombination coefficient, 

/, — = LM(iV,S)- (16) 

Let us point out that any ratio aijaj (particularly, the case 
j - \ and thus ^) is independent of the temperature and de- 
pends exclusively on atomic quantities, namely cross-section, 
transition frequency and statistical weight, a fact which follows 
from Has. lilandim 

v2 



a; 



[rij 



ai_ 



— exp 



h{vi - Vj) 



— . (17) 

J X'j/ 8 J 

Collecting terms, our approximate ionization equilibrium fi- 
nally reads 

' hvi I 1 1 ' 



X a^+CN + CM + Cs) 

= (-) wJ^^I+Cm + Cn + Cs), 



(18) 



where the second variant uses the LTE ratio evaluated at ac- 
tual electron-density and radiation temperature of the ionizing 
continuum. 

4.3. Excitation 

The remaining step concerns the term in the bracket above, i.e., 
the approximate calculation of the excitation inside the lower 
ion (which, of course, is also required in order to calculate the 
partition functions). For consistency, frequencies (energies) are 
still defined with respect to the ionization threshold, i.e., line 
frequencies have to be calculated from v/j = v, -vj > Q instead 
of the usual definition (upper - lower) which would refer to 
excitation energies. 

4.3.1. Meta-stable levels 

We begin with the occupation numbers of meta-stable levels 
which c an be populated via exc ited levels or via the continuum 
(see also lAbbott & Lucvll98^ . 

Population via excited levels. Denoting the excited level by j, 
considering the fact that this excited level is fed by the ground 
state (otherwise it would not exist in our level hierarchy) and 
neglecting collisional processes, ffie population can be approx- 
imated by 



«(^k, 



m e M (;■ € N > m) 



(19) 



(see also Eq. with 6 » 0), where the dilution factor can- 
cels out. In the following, we have to distinguish between two 
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vij ~ v„,j, T,j 



Vmj <V\j ~ Vim, T^X] ~ Ti.ln 



Fig. 2. Population of meta-stable levels via excited ones (see 
text). 



cases: the meta-stable level lies either close to the ground-state 
or close to the excited level, cf. Fig.|2l 

Case a: low lying meta-stable level. The transition frequencies 
of both transitions are fairly equal, vij ^ v,„j, i.e, Ti-jj ^ Tr mj, 
and we find 



ni gi \ kTr,ij) 



(20) 



r.lj 



Note that the population is controlled by the radiation field T^ ij, 
i.e., from frequencies much larger than the "excitation energy". 

Vim- 

Case b: high lying meta-stable level. Now we have v„,j < v\j x 

Vl„„ 



«m gm I rivij 

— « — exp - 

"1 ^1 



__exp,-^ =M .(21) 



and the population depends on the radiation field at (or close 
to) the excitation energy. 

Population via continuum. The third case comprises a popu- 
lation via the continuum which will only be treated in a crude 
approximation, where a correct evaluation will be given later 
(Sect. 14. 5> . If we neglect for the moment the influence of any 
meta-stable and excited levels, we find from Eq. (I18> with 
^ ^ 1, Cn(m,s ) — * 



gm T^r.l 

= — — — exp 



_h 






k 




^ r,m / 



Note that all three cases converge to the identical result 

nm t nm \ r ^ 

> — for Trad = const 

which is quoted bv lAbbott & Lucvl (Il985l) . 



(22) 



(23) 



In order to continue our calculation of Cm, we find from 
Eqs.lHandini 



^1,1X1 a„, gi 



RiK ai g,„ 



I-- 






k 




-^)1 



(24) 



Multiplying by we find that for the three cases Eg s. 1201 

EJandEI 



n\R 











' Vm 


Vl 




Vlm \ 






exp 


k 


\ -^r,m 


Tr,l 


-H 




a,,, 




exp 




^ Vm 


Vl 


-H 


Vim 


a\ 




k 


K Tr,m 


Tr.l 


Ti;lm 



an. 
ai 



(25) 



respectively. As mentioned before, the result for the third case 
(population over continuum) is only a crude approximation, 
which is also evident from the fact that it depends only on 
atomic quantities but not on any radiation temperature. 

4.3.2. Subordinate levels 

Due to our definition of subordinate levels their population can 
be approximated by a two-level-atom Ansatz (between ground- 
state j - I and subordinate level i e N or between meta-stable 
level j e M and subordinate level i e S), such that the popula- 
tion can be expressed by 



ieN(S), ;6 1(M)(26) 



where 5 is the parameter expressing the competition between 
thermaUzation ((J — > 1) and local escape (in Sobolev approxi- 
mation). 



(27) 



e{\-l5)+li 

e is the usual LTE parameter in a two-level atom, 
Cu 



e - 



(28) 



with coUisional de-excitation rate Cji and Einstein-coefficient 
Aji. p is the local escape probability in Sobolev approximation. 



1 TsijJ-) 

and the illuminating radiation field is approximated by 



1 fi 



1 - e"'^^'''' 



(29) 



(30) 



'II. Ts(fi) 

Note, that our approximation i26\ neglects any coupling to the 
continuum inside the resonance zone. By means of Eq. (I17> . 
the individual terms comprising Cn can be calculated from 



— — exp 

a I T,i 



W(l -5i,)exp 



Jl 


f V,- 




k 







(31) 
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whereas the components of Cs are described by 



and 



[rii RuJi^s 



/meM 



ni_Rii. 
n„, R\ 



(32) 



X exp 



hv„. 



W(l-6„,i)expl- 



kTr 



+ 5,„,exp(- — 



(with (M,„/ni) taken from Eas. 1201 1211 and l22l respectively). 
Obviously, the population of subordinate levels is controlled by 
at least three different radiation temperatures (ionization from 
the considered level, ionization from the connected lower level 
and excitation due to line processes). 

4.4. Limiting cases 

In the following, we will consider some limiting cases which 
have to be reproduced by our approach. 

Constant radiation temperature, no collisional excitation 
are the assu mptions underlyin g the description bv lSpringmannI 
1I1997I) and IPuIs et all ( I2OOO ) on the basis of Lucy's unpub- 
lished notes. With Trad = const, meta-stable levels are popu- 
lated via 



- — (V,„ - Vl -H Vl,„ = 0!), 

«1 ^lW,„eM Q'l 



(33) 



independent of the actual feeding mechanism. With 6-0 (only 
radiative line processes), we thus obtain for the population of 
subordinate levels (both / e and i € S\) 



(Hi Riii \ a, 
-L-ii = W—. 



(34) 



Thus, for constant radiation temperatures, it does not play any 
role how the meta-stable levels are populated, and whether 
subordinate levels are connected to the ground-state or to a 
meta-stable level. Only the corresponding recombination co- 
efficient is of importance and the fact that subordinate levels 
suffer from dilution (since they are fed by a diluted radiation 
fiel d), whereas for meta-s table levels this quantity cancels out 
(cf. lAbbott & LucvIITqSSI) . In total, our simplified ionization 
equilibrium then becomes 



rii 



'eN.S 



(35) 



If we define 77 as the fraction of recombination coefficients for 
all meta-stable levels. 



we find 



eM : 



1 



(36) 



(37) 



(38) 



and the ionization equilibrium can be described in a very com- 
pact way. 



n\ \ni /j.^^^ V Trad ^ ' 



(39) 



which indeed is the result of the previous investigations men- 
tioned above. If we further prohibit all ionizations from meta- 
stable and subordinate levels, i.e. allow for 

lonization/recombination only from and to tine ground-state, 
we find with (, - I and 77 = 



T^rad 



which is a well-known result and also valid for the cas e where 
all lines are optically thick and in detailed balance, e.g.. lAbbot3 
a982) . The 



LTE-case is recovered independently from the specific values 
of ^ and 77 in the lowermost atmosphere, when the dilution fac- 
tor approaches unity, W - I, and the radiation field becomes 
Planck, Tiad — * Te. In this case, the ionization balance becomes 



rii 



— = — (^ + 77 + (l-/7-0)= — 



and for the excitation we have 
/ 6M,N,S. 



Hi 



(41) 



(42) 



4.5. Different spin systems 

The last problem to overcome is the presence of different spin 
systems, a problem already encountered for Hei . Our approx- 
imation is to consider the different systems as completely de- 
coupled (except if strong inter-combination lines are present, 
see below), since a coupling via collisional inter-combination 
is effective only at high densities (i.e., in or close to LTE, which 
is treated explicitely in our procedure anyway). 

Then for each of the separate multiplets, the ionization 
equation can be calculated independently. The different sub- 
systems are defined in the following way 

• the first subsystem includes all levels coupled to the 
ground-state plus those meta-stable levels fed from higher 
lying (subordinate) levels (case a/b in Sect. 14.3. 1> . In this 
way, we include also systems of different spin which are 
connected to the ground-state system via strong inter- 
combination lines, a condition which is rarely met. In total, 
the ground-state subsystem includes the levels ; € 1, N, M', 
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S', where M' comprises all case a/b meta-stable levels and 
S' those excited levels which are coupled to M'. For rea- 
sons of convenience, we will denote this set of levels by (1, 
N'). 

• a second group of j subsystems comprises 

- systems of different spin decoupled from the ground- 
state; 

- "normal" meta-stable levels populated via continuum 
processes (poorly approximated so far) and excited lev- 
els coupled to those. 

Both groups can be treated in a similar way and are also 
identified in a similar manner, namely from the condition 
that the lowest state of these systems is meta-stable and not 
fed from higher lying levels. Each subsystem comprises the 
"effective" ground state nij e M" (either different spin or 
fed by continuum) and coupled levels, / e S"j. 

Once more, j is the number of meta-stable levels per ion which 
are not fed by higher lying levels. Note that for a single spin- 
system with meta-stable levels, there are now \ + j different 
subsystems if /' continuum fed meta-stable levels were present. 
Note also that by using this approach we neglect a possible 
coupling of two or more non-ground-state multiplets via strong 
inter-combination lines, if there were any. 

Because of the assumed decoupling, for each subsystem 
we can write down the appropriate ionization equation. For the 
ground-state system, we have 



n„ I nA i Te 

n, Ri, 



Cn' 



■sr-t n i K 



(43) 
(44) 



1 t^u 



where, again, N' comprises the "old" levels € N, M' and S'. 

Note the difference between Zfi and ^ from Eg. 1161 

For each of the j additional subsystems, we obtain in analogy 



fill in,. 



Cnij — 



Ii(m„Sp ffi 



(45) 
(46) 



and S"j comprises all levels coupled to my. The individual com- 
ponents of Cap and Cs" are calculated as described in Sect. 14. 3 
Dividing Eq.|^by Eq.|^ we find for the ratios {n^ ln\) (re- 
quired, e.g., for calculating the partition functions). 



the ground-state plus a number of meta-stable levels alone, we 
would have Cat' = Cs" = 0, f i = 1 and ^,„, = a\la„i.. In 
this case. Eg. 148 land Eg. 1221 would give identical results, which 
shows that both approaches are consistent under the discussed 
conditions. But as already pointed out, Eq.|52is highly approx- 
imative if a variety of levels are involved, and the occupation 
numbers should be calculated according to Eg.l48lalwavs. 

The major difference to our former approach (one spin sys- 
tem only) is the following. In approach "one", the ground-state 
population, n,f/«i, is affected by all meta-stable levels, whereas 
in approach "two" only those meta-stable levels have an influ- 
ence which are coupled to the ground-state system via higher 
levels. 



Constant radiation temperature, no collisional excitation. 
Concerning the limiting case where Trad = const and <5 = 0, 
Eq. remains valid if we account for the different "normal- 
ization", i.e., if we replace f by and include only those meta- 
stable levels into r\ which are populated via excited levels: 



n\ 

with 



^rad 



'71 



(^l+77l+W(l-7?i-^l)) 



(49) 



(50) 



Inside the individual sub-systems we then obtain 



n^ 



w 



- rad 



U'+w(i-nl r = — ^«,/(5i) 



which immediately indicates the correct thermalization for 
W = 1 and Trad ^ T^. After dividing Eq.g^lbylsT] we find 
for the population of (n„./ni) in the same limit 



- ^Xp 

"1 8l \ kTrid 



(' + w(i - n 



(52) 



This expression reveals two things. First, we obtain the correct 
population in LTE when W — > 1. Second, the difference to 
our crude approximation in Sect. 14. 3.T1 becomes obvious: the 
quasi-LTE ratio (I23> has to be multiplied by the last factor in 
the above equation to obtain consistent populations. This factor 
(which can be lower or larger than unity) becomes unity only 
when W — > 1 (i.e., in the lower atmosphere) or for - ^' - 
1, i.e., in those cases where only ground-state and meta-stable 
levels are present, as already discussed above. 



(nJniY 

"ntj ir,l 



"1 iniilnm)* 

r,mj 

or, explicitly written 



exp 



l + Cs 



(47) 



I+Cn' 

l+Cy;j 



.(48) 



The last equation is the "correct approximation" for continuum 
fed meta-stable levels. On the one hand, if the ion consists of 



4.6. Accelerated Lambda Iteration 

To overcome the well-known problems of the Lambda-iteration 
when coupling the rate-equations with the equation of radia- 
tive transfer, we apply the concep t of the Acceler ated Lambda 

I 1 1 1 

Iteration (ALI, for a review see Hubenv 1992) to obtain a 
fast and reliable convergence of the solution. Since our rate- 
equations have been formulated in a non-conventional way 
and since the radiation field is expressed in terms of local, 
frequency-dependent radiation temperatures, the procedure has 
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to be modified somewhat, and we will describe the required re- 
forn iulations as follows (for a comparable implementation see 
also lde Koter et al.ll993l) . 

At first, assume that only one bound-free opacity is present, 
i.e., the radiation-field is controlled by the opacity of the con- 
sidered transition / (no overlapping continua present). In this 
case, the usual ALl formulation for the mean intensity J" at 
iteration cycle n is given by 



= ^Jy + AIS'J with AJy = - ^V^ 



(53) 



where S , is the continuum source-function for transition / and 
A* the corresponding Approximate Lambda Operator (ALO), 
calculated in parallel with the solution of the continuum trans- 
fer^ fo llowing the method suggested by .Rvbicki & Hummer. 
( Il99ll Appendix A). 

Substituting this expression into the rate equations, we find 
for the corresponding effective ionization/recombination rate 
coefficients 



Anuv 
hv 



■AJy dv 



r°° \nay (2hv^ 



(54) 



i.e., the problematic, optically thick part of the radiation field 
has been canceled analytically. Neglecting again stimulated 
emission (the A7v-term in the recombination rate coefficient 
above), approximating 5""' = By/b"^^ with Planck-function 
By and NLTE-departure coefficient ', and using the radia- 
tion temperature at the threshold, Trad.i along with Seaton pa- 
rameters yS = 1, s = 2, we have in analogy to Eg. 1141 



8n , I kT 



n-\ 



A* fcT^e ^_/,y.ytr„ 
b"-^ h 



IK 

——aiVje, (I - Ay) .- R^i - Ri^i 

n 



(56) 
(57) 



In those cases where an overlapping continuum is present, i.e., 
if different transitions contribute to the opacity, the ALO has to 
be modified according to 



Al^PMAl with Pi 



A'tot(v) ' 



(58) 



Xi is the opacity of the considered transition, tot the total opac- 
ity and Pi is assumed to be constant between two subsequent 
iterations (cf. Paper I). Let us mention explicitely that the opac- 
ities used for the radiative transfer are calculated from their 
actual Seaton parameters (/?, s), whereas the uniform values 
ip - \, s - 2) are applied "only" to evaluate the approximate 
ionization/recombination rates. 

Since the Lambda Iteration fails only in the optically thick 
case, we apply the ALI-scheme exclusively for ground state 



* including the pseudo-continua from the multitude of overlapping 
lines, cf. Sect.|5| 



transitions. Thus, by substituting the effective rate coefficients 
R\ii and Ri^i into Eqs.0]|5] we have 



rii 



R\ 



1 _ 111! 4- y Hi-Em. 



n\j Rk\ ■[ + y ("•X ' 



(59) 



Using, again, Eq.^jand the definitions given in J44> . we finally 
obtain 



R^ 



1 



1 



R\. 



\+Cn' Ru 



n\ \ni) R^i 1 _ A fk 

^A-1 



(60) 



In the case of A* s (implying R'^^ = R'^^ - 0), we immedi- 
ately recover the original result, Eq.|^ since 



rill RkI 




(61) 



by means of Eq.[2] If, on the other hand, the ALO is significant 
(i.e., close to unity), we find 



Ik 
Rk\ 



Al 



Wb"-^ T,i"-^ 



— exp 



hvi I 1 



1 



(62) 
(63) 



Thus, the reformulated ALI-scheme collapses to a simple cor- 
rection of the original equation (I43> for the ground-state popu- 
lation. 



— = —(A! = 0) ■ CA(r,- 1"-', b"-\ with factor (64) 



a: 



Ca 



(1 + c^Oww-i 



hvi I 1 



1-^1 a; 



The consistency of this scheme is easily proven, because after 
convergence we would get (cf. Ea.l43> 



1 



(a 








n\ 





W — — exp 



hvi I 1 



^1 (1 + Cn'), 



(65) 



so that the "ALO-correction factor" Ca becomes unity. 
Throughout the iteration the correction factor can take values 
smaller or larger than unity, leading to a fast and reliable con- 
vergence. 

4.7. Test calculations 

In order to check the accuracy of our approximate approach, we 
will present two different test calculations. The first test aims at 
a clean investigation of the methods outlined above, unaffected 
by additional complications such as line-blocking/blanketing. 
To this end, we have computed a pure H/He atmosphere at 
Tss =40,000 K, for two different sets of parameters: the first 
model (A4045 with log g =4.5) corresponds to a dwarf with 
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tau_Ross tau_Ross 

Fig. 3. Approximate NLTE vs. the exact case: He ionization fractions (from top to bottom: Hem, Hen, Hei) for pure H/He 
atmospheric models at T^ff =40,000 K (left panel: dwarf with log §=4.5 and thin wind; right panel: supergiant with logg=3.7 
and thick wind). Bold: exact solution for helium; dotted: He in approximate NLTE (see text). 



thin wind, the second (F4037 with log §=3.7) to a supergiant 
with thick wind.^ 

For both models we have calculated an "exact" solution as 
described in Paper 1, namely by solving for the H/He occupa- 
tion numbers from the complete rate equations, with all lines 
in the CMP and a temperature stratification calculated from 
NLTE Hopf -parameters. In order to test our approach, we cal- 
culated two additional models, with an exact solution for hy- 
drogen only, whereas helium has been treated by means of our 
approximate approach. (In the standard version of our code, 
helium is always treated exactly.) 

Fig EJshows the very good agreement of the resulting ion- 
ization fraction for helium in both cases. The small differences 
at large optical depths (i.e., for LTE conditions) are due to the 
different atomic models for helium used in both the exact and 
the approximate solution. (The data-base applied to the approx- 
imate solution comprises a lower number of levels for both Hei 
and Hen , so that the partition functions are somewhat smaller 
than in the exact case, and consequently also the ionization 
fractions. The occupation numbers of the levels in common are 
identical though). 

Most intriguing is the excellent agreement of the Hen 
ground state departure coefficient as a function of depth (Fig.|4] 
upper panel). The crucial feature is the depopulation of the 
Hen ground-state close to the sonic point, which is a sophis- 
ticate NLTE-effect arising in unified model atmospheres and 
depends on a delicate balance between the conditions at the 
Hen ground-state, the n - 2 ionization edge and the HenLya- 
line (which in itself depends on the radiation field at 303 A 
and the escape probabilities), cf. Gabler et al. ( 1989). The com- 
parison between exact and approximate solution shows clearly 
that our approach, accounting for frequency dependent radia- 



tion temperatures and important line transitions, is actually able 
to cope with such complicated problems.** 

In the lower panel of the figure, we have displayed 
the "ground"-state departure coefficients of Hei, for the 
triplet and singlet system (upper and lower curves, respec- 
tively). Although the precision is not as excellent as for the 
Hen ground-state, one has to consider that Hei at 40,000 K is 
an extremely rare ion, and that the major features (depopula- 
tion of the singlet ground-state, no depopulation for the triplet 
ground-state) are reproduced fairly well. 

The second test investigates the behaviour of the metals. 
We compare the results from the approximate method with re- 
sults from an "almost" exact solution, for model F4037. As we 
will see in Sect. the introduction of a consistent tempera- 
ture structure calculated in parallel with the solution of the rate 
equations forced us to consider the most important elements 
(in terms of their abundance) in a more precise way than de- 
scribed so far, at least if we calculate the temperature from the 
electron thermal balance. In this case it is extremely important 
that the occupation numbers from all excited levels are known 
to a high precision in order to account for the cooling/heating 
by bound-bound collisions in a concise way. Unfortunately, this 
latter constraint cannot be fulfilled by our approximate method, 
simply because not all excited levels are considered, and small 
deviations from the exact solution (which are negligible for the 
effects of line-blocking, see below) can have disastrous effects 
on the total cooling/heating rates. 

Thus, for the most abundant elements the complete set of 
rate-equations has to be solved for in any case, and this solution 
(which uses a Sobolev line transfer, cf. Sect. is compared 
with our approximate one in Fig.|5] for the ionization stages m 
to v of some important metals, namely C, O, Si, Ar, Fe and Ni. 
Note that the comparison includes the effects of line-blocking 
on the radiation field, where this radiation field has been cal- 



Conceming the nomenclature of our models, cf. Sect.|5| 



* Actually, it was this feature that motivated us to refrain from fre- 
quency independent radiation temperatures, since a first comparison 
using the latter simplification gave extremely unsatisfactory results. 
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tau_Ross tau_Ross 



Fig. 5. Approximate NLTE (grey) vs. the results of a solution of the complete rate equations, using Sobolev line transfer (black): 
ionization fractions of important metals for model F4037. Displayed are the ionization stages m, rv, v (dotted, dashed and dashed- 
dotted, respectively). 
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Fig. 6. As Fig.|51 but for model A4045. 
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tau_Ross 

Fig. 4. Approximate NLTE (dotted) vs. the exact case (bold): 
He departure coefficients for model F4037. Upper panel: 
Hen ground-state departure coefficient. Lower panel: Hei triplet 
and singlet "ground"-states (upper and lower curves, respec- 
tively). 

culated either from the exact occupation numbers or from the 
corresponding approximate values, respectively. Our compari- 
son demonstrates three important points. 

- The transition between LTE and NLTE (taking place at 
tr > 2/3 in our approximate approach) is described cor- 
rectly. 

- The approximate treatment works particularly well for el- 
ements with complex electronic structure (Ar, Fe, Ni), i.e., 
our treatment of meta-stable levels is reasonable. 

- If there are differences, they occur predominantly in the 
outer wind. 

In almost aU considered cases, the principal run of the approx- 
imate ionization fractions agrees reasonably or even perfectly 
well with the exact result. The only exception is oxygen where 
the major/minor stages (iv/iii) appear reversed in the outer wind 
(no problems have been found for nitrogen and neon which 
are not displayed here). These differences in the outer wind 
(see also C and Si) are partly due to two effects. On the one 
hand, our approach becomes questionable in those cases when 
all hne transitions are optically thin so that the two-level-atom 




1.5x10 t , , , , i 

100 1000 10000 

ambda [A] 

Fig. 7. As Fig.lS] Comparison of radiation temperatures of con- 
verged models. 

approach fails to describe the excitation-balance of subordinate 
levels. If only this effect were responsible this would imply (as 
suggested by our referee) that the discrepancy should become 
worse for thinner winds. Thus, we performed a similar compar- 
ison for model A4045, which has a considerably lower wind 
density than model F4037, by a factor of almost 100. The cor- 
responding ionization fractions are shown in Fig. |6| Note that 
the transition point between photosphere and wind is located 
at lower values of tr, compared to model F4037, due to the 
weaker wind. Interestingly, the discrepancies between approx- 
imated and "exact" ionization fractions in the outer wind have 
remained at the same level as for model F4037, and in the case 
of oxygen the situation is almost perfect now. Consequently, 
the effect discussed above cannot be responsible alone for the 
observed discrepancy, and we attribute it to a combination of 
various ingredients inherent to our approximative approach. 

For our models, however, this is of minor importance, since 
we are not aiming at a perfect description of the occupation 
numbers in the outer wind unless we actually need it, i.e., when 
a consistent temperature structure shall be derived. In this latter 
case, the occupation numbers are calculated exactly anyway. 

Different occupation numbers influence the radiation field, 
which in turn influences the occupation numbers, and so on. 
This is the second important process which might affect our 
final approximate solution. Fig.0compares the emergent fluxes 
(expressed as radiation temperatures) for the converged models 
of F4037, calculated by both alternative approaches. 

Due to the excellent agreement between the ionization frac- 
tions in the line/continuum forming part of the atmosphere, also 
the fluxes agree very well. The maximum differences, located 
between 200 to 400 A, are of the order of +1,000 K, which 
translates to a typical difference in population of ±0.15 dex in 
the outer wind. 

Globally, however, the differences in flux are so small that 
we can consider the two results as equivalent. Thus, the radi- 
ation field calculated in parallel with the line-blocking back- 
ground elements is insensitive to the chosen approach (exact 
vs. approximate occupation numbers) which primarily differs 
in the precision (and presence) of subordinate levels. 
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5. Approximate line-blocking 

The most time-consuming part in the computation of realistic 
stellar atmospheres is the calculation of the radiation field, re- 
alizing the multitude of overlappi ng^ lines with co nsiderable 
opacity (see also the discussion bv lPuls & Pauldrachfcl99Q. and 
[Pauldrach et al. 2001). 

For CMFGEN as well as for the wind-code de veloped by the 
Potsd am group (for a recent status report, see lOrafener et al.l 
l2002h . this problem has been tackled by performing a 
comoving-frame solution for the complete EUV/UV range. 
Obviously, this approach is very time-consuming. A quick cal- 
culation shows that the number of frequency points which must 
be treated is of the order of 900,000, if a range between 200 and 
2,000 A and a typical resolution of 0.8 km s^'is considered 
(i.e., ten points covering a thermal width of 8 km s '). 

In the approach followed by WM-Basic, on the other hand, 
an observer's frame solution is performed which requires 
"only" a couple of thousand frequency points to be considered. 
The conservation of work, however, immediately implies that 
in this case a lot of time has to be spent on the resolution of the 
resonance zones of the overlapping lines, a problem which is 
avoided a priori in a CMF calculation. 

In order to solve the problem on a minimum time-scale, 
both a Monte-Carlo solution"* (e.g., [Schaerer & Schmutz 
lll994 : ISchaerer & de KoteH lll997l) , and a statistical approach 
are feasible. 

Since the number of metal lines to be treated is very large, 
the information about the exact position of individual lines in- 
side a (continuum transfer) frequency grid interval becomes 
less importan t for obtaining a repr esentative mean background. 
As shown bv IWehrse et alJ lll998 ). the Poisson Point Process 
is well suited to describe such a line ensemble, particularly be- 
cause it is very flexible and can be described by relatively few 
parameters. 

The additional introductio n of a Generalized Opacity 
Distribution Function by Basche lTet al.l dioOl) serves two pur- 
poses. First, additional analytical insight is given into the ef- 
fects of the vast amount of blocking lines on the mean opacity 
in differentially moving media with line overlap. Second, it is 
a fast tool to derive such mean backgrounds numerically. In 
particular, it is able to "solve problems that have been inacces- 
sible up to now as e^. the influence of very many, very weak 
lines' 'te aschek etalfcOOll) . and to describe the transition from 
a static to a moving configuration, since it is equally efficient 
in both cases. 

In our opinion, this approach is very promising, and work 
adapting and applying the corresponding method is presently 
under way in our group. Since it will take some time to finalize 
this approach (the most cumbersome problem is the formula- 
tion of consistent emissivities), we have followed a somewhat 
simplified approach in the mean time, which relies on similar 
arguments and has been developed by carefully comparing with 
results from "exact" methods, mostly with the model g rid cal- 
culated with WM-Basic as described bv lPauldrach et al T iEoOl) . 



Again, the principal idea is to define suitably averaged 
quantities which represent a mean background and which can 
be calculated easily and fast. The multitude of lines will be 
approximated in terms of a pseudo continuum (split into a 
"true" absorption and a scattering component), so that the ra- 
diative transfer can be performed by means of a standard con- 
tinuum solution, for relatively few frequency points (see be- 
low). Strongest emphasis has been given to the requirement that 
any integral quantity calculated from the radiation field (such as 
the photo-integrals) has to give good approximations compared 
to the exact case, because these quantities (and not the frequen- 
tial ones) are most decisive for a correct description of the level 
populations and, in turn, for the blocked radiation field. 

5.1. Mean opacities 

To this end, we define a "coarse grid" with spacing 2A^v'™, 
where vJH is a typical thermal velocity (say, of oxygen) in- 
cluding micro-turbulence, and 2N is an integer of the order of 
100. (The reason to define here 2N instead of will soon be- 
come clear.) Under typical conditions, this grid has a resolu- 
tion of 1,000-1,500 km s^'and is used to calculate appropri- 
ate averaged opacities. With respect to a simplified approach, a 
mean constructed in analogy to the Rosseland mean is perfectly 
suited, i.e., an average of the inverse of the opacity. 



f 

1 JlNv'i 



dv 

.tot 



/ 



(66) 



dv 



' both in the observer's and in the comoving frame 
which becomes costly as well if a detailed description of all pos- 
sible interactions between radiation field and plasma is accounted for. 



since it has the following advantageous properties: 

a) if no lines are present, the pure continuum opacity is recov- 
ered 

b) if one frequency interval is completely filled with non- 
overlapping, strong lines of equal strength, also the average 
opacity approaches this value, whereas 

c) in those cases when the interval has "gaps" in the opacity, 
these gaps lead to a significant reduction of the mean, i.e., 
allow for an appropriate escape of photons. Note that any 
linear average has the effect that one strong line alone (of 
typical width 2vJJj) would give rise to a rather large mean 
opacity (just a factor of weaker than in case b) and, thus, 
would forbid the actual escape. 

d) Finally, the average according to Eq.|5^is consistent with 
the standard Rosseland mean in the lowermost atmosphere 
(as long as dBy/dT is roughly constant over one interval), 
i.e., it is consistent with the diffusion approximation applied 
as a lower boundary condition in the equation of (contin- 
uum) radiative transfer. 

Because of the large number of contributing lines (typically 
5-10^ (O-type) to 10^ (A-type) lines if only ions of significant 
population are considered"), the calculation of this mean has 
to be fast. 



" Remember, that our present data base comprises 4.2 • 10^ lines in 
total. 
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First, assume that any velocity field effects (leading to 
Doppler-shift induced line overlaps) are insignificant, i.e., as- 
sume a thin wind, so that line blocking is essential only in the 
subsonic regions of the wind. The generalization in order to ap- 
proximate line-overlap in the wind will be described later on. 

Instead of evaluating the "exact" profile function, for each 
line we use a box car profile of width 2v^. The frequential line 
opacity is, thus, given by 



XL 



Xh for vo + Avd < V < vo - Avd 
else 



1 ne rii 

g/-, 

2Avd nieC gi 



Avn = 



(67) 
(68) 



where stimulated emission has been neglected again. Due to 
this definition, at least the frequency-integrated line opacity is 
correctly recovered. The coarse frequency grid is now divided 
into sub-intervals of width Av - 2Avd. Inside each of these 
sub-intervals ("channels") we sum up any line opacity which 
has appropriate rest-wavelength. Insofar, we account (approxi- 
matively) for any intrinsic (i.e., not wind-induced) line overlap. 
Inside each channel /, we thus have a (total) frequential opacity 



■i - Z^Xlj +Xv 



xT-xT' + <T. 



(69) 



if lines j are located inside channel / and the continuum opac- 
ity is assumed to be constant inside each coarse grid interval. 
X^v^'^ is the contribution by true absorption processes, and a-f. 
the contribution by electron scattering. After replacing the inte- 
grals by appropriate sums and since all channels have the same 
width, the mean opacity (on the coarse grid) is simply given by 



^ A 
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2^ Ay 
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(70) 



For later purposes we split this mean opacity into the contri- 
bution from lines and continuum, respectively, where the line- 
contribution is given by 



(Xl)^ 



N 



Xv 



(71) 



and we have 



(X'°')-(XL}+Xy 



(72) 



Note that both mean opacities, (xl) and are frequency 

dependent as a function of coarse grid index. In accordance 
with our reasoning from above, Eq. (17 1> implies that 

a) ifxLj - for all lines inside one interval, the correct result 
(Yl) = is obtained 

b) if the same total line opacity Xhiv) is present inside each 
channel, this value will also be obtained for the mean, 
CTl) =Xl(v). 



c) if only one (strong) line is present, the mean line opac- 
ity is given by (xi) ^ ^™°'/(A^ - 1). i-S-^ it will be much 
smaller than the continuum opacity, since most of the flux 
can escape via the (A^ - 1) unblocked channels (according 
to our present assumption that Doppler-induced line over- 
lap is negligible). 

Finally, let us point out that the opacities constructed in this 
way are used also to calculate the photospheric line pressure, 
in analogy to the description given in Paper I (Eq. 3), however 
including the line contribution (cf. Fig. II U .'^ 

5.2. Emissivities 

In order to calculate the corresponding emissivities, we assume 
that each transition can be described by means of a two-level 
atom, where the lower occupation number is known from the 
solution ("exact" or approximate) of the rate equations.'^ 

Although this assumption is hardly justified for (weak) re- 
combination lines, it is a fair representation for most of the 
stronger transitions arising from either the ground-state or a 
meta-stable level, particularly if the level population itself is 
calculated from a multi-level atom. 

It might be argued that the two-level atom approach is su- 
perfluous for those connecting transitions which are calculated 
from an exact NLTE solution, since the occupation numbers for 
both levels and, thus, the source-functions are already known. 
The maximum number of these lines is of the order of 30,000, 
and therefore much lower than the total number of lines we 
are using for our Une-blocking calculations (cf. Sect. O. For 
the latter transitions, however, only the lower level is present in 
the atomic models, so that the corresponding source-functions 
have to be approximated in any case. 

Moreover, treating all lines (including the connecting tran- 
sitions) in a two-level way has the additional advantage that 
the contribution of scattering and thermal processes can be eas- 
ily split, which allows to simulate their impact by means of a 
pseudo-continuum, so that the standard continuum transfer can 
be applied without any modification. 

To keep things simple and as fast as possible and to be in 
accordance with our assumption of box car profiles, we replace 
the scattering integral inside the two-level source-function by 
mean intensities, i.e., we write 



Sh - pJv + SBv, p — 1-6, 



(73) 



where 6 has been defined in Eq.|^and is evaluated for the line- 
specific thermalization parameter and escape probability. The 



In our present version of fastwind we allow for deviations from 
the generalized Kramer-law (Paper I, Eq. 2) by simply including 
theses deviations as correction-factors into the atmospheric structure 
equations. This method becomes important for models at rather cool 
temperatures when hydrogen and background-metals are recombining 
(and become ionized again) in photospheric regions, which usually 
leads to some deviations from the above (power-) law. 

Note that this approach is equivalent to the typical assumption 
made if deriving the radiation field via Monte-Carlo simulations. 
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total source-function (in channel /, before averaging) is then 
given by 



S v,i — 



T/^'""" + O-e/v + IjjXhj (PjJv +SjBy) 

. .cont _i_ V ' 
Xv + LjXhj 



(74) 



with 77v'™'' being the thermal component of the continuum 
emissivity. Note that the frequential line-opacity xlj includes 
the "profile function" (2Avd)"', cf. Eg. 1671 

In the following, we will investigate how to average the 
above quantities in order to be consistent with our definition of 
(^'°') and (xl)- With respect to the equation of transfer, which 
will be finally solved on the coarse grid, we find that after inte- 
gration over the subgrid-channels 



1 



(75) 



with z being the depth variable along the impact parameter p in 
the usual {p, z)-geometry. Strictly speaking, the first term in the 
above equation (i.e., the mean inverse opacity) is given by 



1 fHM^ 



(76) 



(where the denominator is equivalent to d/dz(/v), and all inte- 
grals extend over the range 2A^vJJj), i.e., a different definition ap- 
phes when compared to the corresponding quantity in Eq. l66l 
Our crucial approximation is to equate both definitions, i.e., in- 
side each coarse grid cell (of width » 1,000. . . 1,500 km s"') 
we assume that 



rid r d r dv r 

--lydv/ -lydV ^ — / 

J Xvdz J az J Xv J 



dv. 



Let us frankly admit that this approximation can be justified 
only if a) the spatial gradient of the specific intensity is a slowly 
varying function of frequency (i.e., deep in the atmosphere) or 
b) the opacities are similar for most of the sub-channels, i.e., 
either no lines are present at all or the (summed) line-opacities 
do not vary too much. Additionally and most important, this 
approximation still works in those cases when only a couple 
of channels are populated by large opacities and the rest is 
filled by a weak background due to the inverse relation be- 
tween opacity and intensity: on the Uis, the high opacity chan- 
nels do not contribute to the fraction because of the correspond- 
ingly low intensities in both the nominator and the denomina- 
tor, whereas on the rhs these channels drop out at least in the 
nominator because of the low value of Ijx- 

There are, of course, a number of cases where the above ap- 
proximation is only poor. With respect to the results presented 
below and since we are not aiming at a perfect, highly resolved 
description of the radiation field in the line-blocking EUV/UV 
regime, the errors introduced by the above approximation (and 
the following one, which is of similar quality) are acceptable 
though. 

In order to proceed with appropriate expressions for the 
emissivity, the mean source-function, {S v), is given by 



where the Planck-function By is assumed to be constant within 
one macro-grid interval. For those averages which are multi- 
plied by (Jy), we have employed an approximation similar to 
the one discussed above. If we, finally, denote the opacity de- 
pendent means of the third and fourth term by /nth and /h, re- 
spectively (i.e., non-f/zermal/f/iermal), the equation of radiative 
transfer for the averaged quantities becomes 



+ (cr, + (y'°')/nth) (Jy) + (y'°')/th By 



(/v>, (78) 



and can be solved in the conventional way (pure continuum 
transport). The resulting quantities for the radiation field are to 
be understood as average quantities, in the sense that integral 
quantities such as J Jydv or j Hydv are described correctly, at 
least in most cases. The coefficients f± and /nth can be calcu- 
lated by summing over the sub-channels. 
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(80) 
(81) 



and after some simple algebraic manipulations, the following 
relation is obtained: 



CrL) = (;r'°')(/th + /nth). 



(82) 



Xv 



Xv 



With this equation it is easy to show that the mean source- 
function ( I77> allows for a correct thermalization, if rf/™^ — > 
^c.trae^^ and {Jy) — > By. In this case, the mean source-function 
becomes {Sy) - By, q.e.d. 

We now need to incorporate the effects of the velocity field 
into our approach. Due to the method to average the opac- 
ity, we cannot simply shift the lines with respect to the stel- 
lar frame. Consider, e.g., one strong line to be present without 
any other interfering lines. In "reality" and in the observer's 
frame, the absorption part of this line becomes broader as a 
function of velocity, i.e., the larger the velocity the more flux is 
blocked (of course, a significant part is reemitted due to scat- 
tering). If we simply shift our line(s) as a function of veloc- 
ity, almost nothing would happen, since, as shown above, the 
mean opacity/radiation field remains almost unaffected by one 
strong line, due to the possible escape via the (N-1) unblocked 
sub-channels. Thus, in order to simulate the physical process, 
we proceed in a different way. When the velocity shift becomes 
larger than twice the average "thermal" width (including micro- 
turbulence), we combine (in proportion to the local velocity) 
more and more subchannels to increase the relative weight of 
the line in the mean opacity. In particular, the line width (more 
precisely, the width of the sub-channels) is set to the value 

Av = max(2AvD, — ). (83) 
/lo 

Although this procedure is highly approximative, it allows to 
deal with the effects of "line-shadowing" and prevents any pre- 
mature escape of photons when the lines begin to overlap. 




Fig. 8. FASTWiND VS. WM-Basic (gr ey): comparison of emergent fluxes for two dwarf and two supergiant models at 35 and 45 kK 
(for parameters, cf. IPauldrach et al. 2001) . In order to allow for a meaningful comparison, the high resolution frequency grid 
provided by WM-Basic has been re-mapped while keeping the corresponding flux-integrals conserved. 



5.3. Tests/Comparison with wm -Bas/c 

Before we test our approximate approach by comparing with 
alternative calculations, let us mention two important consis- 
tency checks we have performed. 

a) The calculated models (and spectral energy distribu- 
tions/line profiles) are, almost, independent of the actual 
value of coarse grid cells, A^, at least if varied within a rea- 
sonable range. (We checked for values between 0.5A^ to 2N, 
for N = 60.) 

b) As long as the IR/radio-range is not considered, our sim- 
pler models with a temperature-structure calculated from 
Hopf-parameters and all background elements in approx- 
imate NLTE agree very well with complex models in- 
cluding a consistent T-structure. This check verifies anal- 
yses performed with previous versions of fas twind, e.g., 
lHerreroetai](l2002h and.Repolu st et al.l ( l2004 . 

In the following, we will compare the fluxes from our models 
with those calculated by WM-Basic by means of the O-star grid 



presented bv'Pauldr ach et alJ fcOOlV^. These tests should give 
reasonable agreement, since both codes use the same atomic 
data base for the background-elements. A comparisons with 
results from cmfgen will be discussed later on. 

The parameters of the corresponding models (calc ulated 
without X-rays) can be found in [Pauldrach et all (l200ll 
Table 5). Our models have been constructed as closely as pos- 
sible to the approach inherent to WM-Basic, i.e., including a 
consistent temperature stratification (which will be described 
in Sect. and Sobolev line-transfer. For the velocity-field, we 
have used /3 - 0.9, which results in a stratifiction very close to 
the one predicted by WM-Basic (see below). The computation 
time on a 2 GHz processor machine is of the order of 15 to 20 
minutes per model (typically 40 to 50 iteration cycles for a final 
convergence below 0.003 in all quantities, if the temperature is 
updated each 2nd cycle). 

The grid comprises 6 "dwarfs" and 5 "supergiants" in the 
range between 30. . . 55 kK ("D30". . . "D55" and "S30". . . 
"S50", respectively), and we have concentrated on the grid with 



available via http://www.usm. uni-muenchen.de/people/adi/Models/Mode 
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Fig. 9. Comparison of ionizing photon number for the model grid provided by IPauldrach et all (l200 ll) . Left panel: logarithm 
of Zanstra-integrals, log gi, (Eg. I84>. for Hi, Hei, On and Hen. Bold/dashed: dwarfs/supergiants as calculated by WM-Basic; 
dotted/dashed-dotted: results from fastwind. Right panel: ratio of corresponding Zanstra integrals, Alog^x = log 2*^^ - 
log 2™ (WMB: WM-Basic, FW: fastwind), for dwarfs (triangles) and supergiants. For the supergiant model at 30,000 K ("S30"), 
we have used the fastwind model without photospheric line-pressure in both figures. The corresponding results for the "correct" 
model, i.e., including photospheric line-pressure, are indicated by the "-h" sign (see text). 



solar abundances, in order to deal with more prominent effects 
related to line-blocking/blanketing. Fig |8] compares the emer- 
gent fluxes for some typical cases, two dwarf and two super- 
giant models at 35 and 45 kK. In order to allow for a mean- 
ingful comparison, we have re-mapped the high resolution fre- 
quency grid provided by WM-Basic while keeping the corre- 
sponding flux-integrals conserved. 

Overall, the agreement is rather good; in particular, the 
range above 400 A is reproduced very well, except for some 
strong absorption/re-emission features which are missing in 
our mean-opacity approach. We have convinced ourselves that 
in all cases also the IR-fluxes (not displayed here) agree per- 
fectly, i.e., the IR flux-excess induced by the wind is re- 
produced equally well in both codes. Major diff'erences are 
"only" present in two regions: most models difl'er in fluxes be- 



low 200 A, although the strength of Hen Lyman-jump itself 
is very similar Mostly, this problem is related to the enor- 
mous bound-free opacity provided by Oiv (and Fev or Civ for 
the hotter or cooler objects, respectively) leading to an opti- 
cally thick wind from the outermost radius point on (in our 
case, Rmax - 1207?*), so that the flux is rather badly defined in 
this frequential region. As we will see from a comparison with 
models calculated with cmfgen (Figs.lTst. these models predict 
a third alternative for A < 200 A, and even the Lyman-jump 
is different. As a result, we consider the ionizing fluxes in this 
wavelength range as not particularly reUable. Moreover, the in- 
fluence of X-rays becomes decisive, implying that any tool for 
nebula diagnostics should use these numbers only with care. 

The second inconsistency is found in the region between 
300 to 400 A. Although this range poses no problem for super- 
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Fig. 10. Comparison of velocity/density structure for model S30. Grey: WM-Basic; bold/dashed: fastwind with/without photo- 
spheric line-pressure, respectively. 



giants, the flux-blocking predicted by fastwind /or dwarfs be- 
tween 35 to 45 kK is larger than calculated by WM-Basic, with 
a maximum discrepancy around 35 kK. The reader might note 
that CMFGEN again produces somewhat different results in this 
range: agreement with WM-Basic is found for dwarfs, whereas 
the fluxes emitted from supergiants are larger compared to both 
FASTWIND and WM-Basic. 

This dilemma becomes particularly obvious if we consider 
the corresponding Zanstra-integrals, 

-^dv, (84) 

which are proportional to the emitted number of ionizing pho- 
tons. In the left panel of Fig. |9] we compare the logarithm 
of g I, evaluated for Hi , Hei , On and Hen , whereas the cor- 
responding ratios, A log 2^ = \ogQf^^ - log (WMB: 
WM-Basic, FW: fastwind), are displayed in the right panel. 
Obviously both codes predict the same numbers in the hydro- 
gen Lyman and in the Hei continuum. As already discussed, 
the situation is much less satisfactory for the Hen continuum, 
where the differences are particularly significant for super- 
giants. Note, however, that the principal dependence of Queii 
on spectral type and luminosity class, which shows the largest 
variation throughout the spectrum (lower left panel), is much 
more consistent than one might expect on basis of the right 
panel alone. In the On continuum (A < 352 A), finally, the dif- 
ferences for the dwarfs at intermediate spectral type are evi- 
dent. 

Note that in this wavelength range the line-density is very 
large, and differences in the treatment of the weakest back- 
ground opacities might explain the established disagreement. 
An argument in support of this hypothesis is given by the 
fact that FASTWIND recovers the results by WM-Basic perfectly 
if a line-list is used which has significantly less (overlapping) 
weaker lines in the considered interval. For a final statement, 
however, more tests are certainly required. Note that a com- 
parison with CMFGEN addressing this point will not solve the 
problem, since the number of fines included in this code is 



mostly lower than described here, because cmfgen uses only 
those lines where the occupation numbers of both levels are 
known, in contrast to our approach which uses also lines where 
the upper level is lying too high to be included into the rate 
equations. 

One last point we would like to mention concerns model 
S30. In a first comparison, we immediately encountered the 
problem that particularly this model provided fluxes which 
showed significantly less agreement at all frequencies than the 
other models (indicated by the plus-signs in Fig.|9}. Comparing 
the models themselves, it turned out that temperature, den- 
sity and velocity structure showed a severe mismatch in pho- 
tospheric regions (cf. Fig. [TO] grey vs. black curves). After 
some tests, we found that both models agree well if the photo- 
spheric line pressure is neglected in fastwind (grey vs. dashed 
curves in Fig. I1Q> . Most likely, this problem is related to the 
treatment of the line pressure in WM-Basic. Whereas the con- 
tinuum forces are calculated from correctly evaluated opaci- 
ties, the line pressure, independent of location, is calculated in 
terms of the force-multiplier concept, utilizing the Sobolev ap- 
proximation. Particularly, §|™ oc t^", with "depth parameter" 
t oc p/(dv/dr). Thus, g^.^^ decreases rapidly in photospheric re- 
gions when the density is large and the velocity gradient small. 

In those cases where the (static) line pressure is non- 
negligible in photospheric regions, the chances are high that 
the above approximation leads to a too large effective grav- 
ity, i.e., too high densities. Actually, this problem is already 
kn own for a long t ime an d has been discussed in fair detail 
in iPauldrach et al particularly Fig. 6c). The reason 

that this problem occurs only in S30 results from fact that 
the Eddington factor is considerably higher than for almost all 
other models (F = 0.52). Insofar, the photospheric line pressure 
has much more impact than for models with either high grav- 
ity or low r. Moreover, at an effective temperature of 30 kK, 
Feiv with its enormous number of lines spread throughout the 
spectrum is the dominant (or almost dominant) ionization stage 
in the "middle" photosphere, thus, contributing a much larger 
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Fig. 11. Comparison of total photospheric radiative accelera- 
tion for model S30 (bold) vs. results from an analogous hydro- 
static Kurucz-model (dotted). Note that the gravitational accel- 
eration for this model is 1000 cms"^, i.e., the radiative accel- 
eration is very close to this value and, thus, of extreme impor- 
tance (cf. Fig.llO>. The deviations at largest depths are due to 
the fact that this model becomes (spherically) extended in the 
lowermost photosphere, an effect which cannot be treated in a 
plane-parallel approach (cf. Paper I). 

amount of static line pressure than for hotter temperatures, 
where Fev or even Fevi are contributing. 

Note that we have also compared our (cooler) models (from 
our grids as described in Sect. O and from additional A-star 
models) with corresponding Kurucz models, where in most 
cases a very good agreement regarding the photospheric radia- 
tive acceleration has been found, e.g.. Fig. Only for mod- 
els cooler than 9,000 K a mismatch becomes obvious, where 
"our" radiation pressure is too low, due to a number of missing 
Fell lines in the optical (improvements under way!). 

In order to allow for a meaningful comparison concern- 
ing our approximate line-blocking, in Fig.|9]we have used the 
results from our S30 model without photospheric line-force, 
whereas the results from the "actual" model (including §^™) are 
indicated by "-h". Independently, however, Fig.[lOl(left panel) 
also shows the validity of our treatment of the transition zone 
from photosphere to wind (cf. Paper I), since in this region 
both velocity fields agree perfectly. (Remember that WM-Basic 
solves the hydrodynamical equations in a consistent way). 

6. Treatment of inverted levels 

One of the more complex problems when solving the cou- 
pled equations of statistical equilibrium and radiative trans- 
fer is the presence of population inversions, which often occur 
in the outermost layers of hot expanding stellar atmospheres. 
The amount of the overpopulation (i.e., > rti/gi) is usu- 

ally small, but even in this case it invokes a number of prob- 
lems concerning the solution of the radiative transfer equation. 
Particularly with respect to the usual concept of using source 
functions, a problem occurs in the transition zone between 
"normal" population and overpopulation, where the source 



function formally diverges. In addition, factors like exp(-T) 
may produce numerical problems for r < 0. In a number of 
codes, this problem is "solved" by setting the upper level into 
LTE with respect to the lower one or by other approximations. 
Since level inversions are particularly present between levels 
responsible for IR-lines and since fastwind aims at a reliable 
solution also in these cases, we cannot afford such approxima- 
tions and have to solve the "exact" case which in turn has an 
influence on the degree of overpopulation itself. In this section, 
we briefly describe how we have solved the problem in fast- 
wind both with respect to the Sobolev approach and within the 
CMF-transport. 

6.1. Treatment of inversions in the Sobolev 
approximation 

Since the Sobolev approach uses only local quantities, a di- 
vergence of the source function is not possible, except for the 
extremely unlikely case that upper and lower populations, nor- 
malized to the appropriate statistical weight, are numerically 
identical. Thus, we can retain the standard concept (optical 
depth and source function) and follow the approach described 
in Taresch et al. ( 1997): in case of a level inversion, the interac- 
tion function [/(Ts,y6p)'^ is split into two parts in order to avoid 
numerical problems. 



U^Ui + Ui, with f/i = l-yS. 



(85) 



/3 is the usual escape probability in Sobolev approximation 
(Eq.l29>. which for the case of inversion is given by 



exp |t^ I - 1 
Its I 



(86) 



and U2 has been described in 'Tare sch et all ( I I997 I. Eq. A13). 

For \/3p\ — > 00, i.e., dominating continuum, U2 approaches 
zero. In the case of dominant line processes, on the other hand, 
and Ts < 0, U2 approaches (j3-l) and U goes to zero. Thus, we 
recover the "classical" result by Sobolev, where the influence 
of continua has been neglected. 

In our app roach, w e have significantly extended the grid 
used bv Taresch et alJ il991.) from which U2 is calculated by 
means of interpolation. Due to the diflferent behaviour of this 
function in different regions of the {ts,/3p) plane, (four) differ- 
ent tables with different degrees of resolution have been calcu- 
lated. The boundaries of the complete grid comprise the area 
between -6 < loglr^l < 2.8 and -6 < log|/3p| < 6. Beyond 
the boundaries, U2 is calculated analytically (by either consid- 
ering the appropriate limits or using a first order expansion). In 
particular, it turns out that 



U2 = 0.5|T5| 

U2 

U2 



log |t5 I < -6 independent of ySp (87) 
log|;0p|<-6 independent of Ts (88) 



1 



log ^p| > 6, log Its | > -6 



(89) 



which describes the interaction between line and continuum pro- 
cesses, where ts is the Sobolev optical depth and f5p the ratio of contin- 
uum to line opacity in a frequency interval corresponding to the ther- 
mal Doppler width, cf . iHummer i&Ribic3 il9ij;|Puis&Hu^e3 
>1988l) . 
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and the limits for large log|Ts| > 2.8 and -6 < log|y6p| < 6 
result from appropriate extrapolations from the pre-calculated 
tables. 

6.2. Treatment of inversions in the comoving frame 

In the CMF solution, the problem of source-function diver- 
gence is inevitable when a population inversion occurs and the 
standard formalism is used. Even if the local quantities are not 
diverging, there will be an implicit divergence just between the 
two depth-points before and at the beginning of overpopula- 
tion, which, due to the applied discretization, will not be han- 
dled consistently. To avoid this problem, it is more suitable to 
work directly with emissivities and opacities rather than with 
optical depths and source functions. Thus, in the case of inver- 
sion, we solve the two coupled equations of radiative transfer 
in the comoving frame according to 



du dv 

dz dx 

dv du 

dz dx 



-XV 



(90) 
(91) 



where u and v are the usual Feautrier variables, x is the fre- 
quency measured from the center of the line in Doppler width 
units, and z is the depth variable along the impact parame- 
ter The opacity is x - Xeiz) + Xl{z, x) and the emissivity is 
J] = J](z, x). 

In order to discretezise the equations with respect to z and x, 
a fully implicit scheme is used. As was shown by Mihalas et al. 
Appendix B) this method is unconditionally stable. 

6.3. Tests 

A number of tests have been performed concerning both the 
Sobolev and the CMF implementations. Most importantly, we 
have also tested models where the above discretization of the 
CMF equations with respect to z has been used for all tran- 
sitions, not only for the "inverted" ones, and found satisfac- 
tory agreement with our standard implementation using a dis- 
cretization with respect to t. 

After convincing ourselves that the algorithms are working 
in principal, we have tested our improved methods by compar- 
ing them with older results (where in case of inversion the up- 
per level and the line source function where set to zero). This 
comparison has been performed for the O-star grid described 
in the previous section. The results were very satisfying, and 
a number of convergence problems originating from the older 
treatment of inverted populations are no longer present. 

The differences in the resulting H/He line profiles (both in 
the optical and in the IR) turned out to be rather small, since 
for our grid parameters these lines are formed below those re- 
gions where the inversion sets in. However, we like to point out 
that a consistent treatment might be important for winds with 
more extreme mass-loss rates and for a number of metallic IR 
transitions with an inversion already occurring in photospheric 
regions. 



a typical example is the Siiv IR transition 4d ^D"" - 4f 



7. Temperature stratification 

As has been previously mentioned, the present version of fast- 
wind allows for the calculation of a consistent temperature 
stratification, utilizing a flux-correction method i n the lower 
wind and the thermal balance of electrons (cf. iKubat et'all 
1999) in the outer part.''' The region where both methods are 
connected is somewhat dependent on mass-loss, but typically 
lies at tr = 0.5. Although the implementation of this method 
is straightforward, and the contributio n of individual pro cesses 
have been discussed in fair detail bv lDrewllll985[ll989l) . three 
points are worth mentioning. 

In order to calculate the appropriate heating/cooling rates 
resulting from collisional bound-bound transitions, the popula- 
tion of excited levels is as important as the population of ground 
and meta-stable ones. This can readily be seen from the fact 
that the net heating rate from a collisional transition between 
lower level / and upper level u can be expressed as 

Qui - Qiu = (n«C„/ - n;C/„) hv,u = niCu,hvu, - 1 j > (92) 

with collisional rates C,,/ and NLTE departure coefficients 
bi, bu- Thus, the ratio of departure coefficients controls whether 
a certain transition heats or cools the plasma and its deviation 
from unity controls the degree of energy transfer Heating re- 
sults from transitions with an upper level being overpopulated 
with respect to the lower one, and cooling vice versa. Thus, the 
occupation numbers of all ionic levels have to be known with 
some precision, and we have to modify our approach when the 
electron thermal balance is used to calculate the temperature 
profile. The approximate NLTE solution as described in Sect.|3] 
simply does not yield the required occupation numbers of ex- 
cited levels (except those which are directly connected to the 
ground or meta-stable level), and any brute force approxima- 
tion would give incorrect heating/cooling rates. 

To overcome this dilemma we incorporated a detailed so- 
lution of the statistical equilibrium at least for those elements 
with large contributions to the net heating rates (positive or 
negative). After some experiments it turned out that the inclu- 
sion of the most abundant background elements C, N, O, Ne, 
Mg, Si, S, Ar, Fe, Ni (plus the explicit elements, of course) is 
sufficient to stabilize the results. For these elements then, the 
complete rate-equations are solved with line transitions treated 
in Sobolev approximation, whereas for the remaining ones the 
approximate NLTE solution is employed. 

The second point to be mentioned regards the flux- 
conservation of the final models. The conventional approach 
to calculate the energy balance, formulated in terms of radia- 
tive equilibrium, satisfies this constraint by construction, at 
least in principle. (Most numerical codes, including cmfgen and 
FASTWiND, calculate mean intensity and flux on different grids, 
which somewhat destroys the coupling between radiative equi- 
librium and flux conservation). On the other hand, our formu- 
lation in terms of the electron thermal balance is decoupled 
from the latter requirement, at least regarding any explicit de- 
pendence. Note, however, that there is an implicit coupling via 

Note that adiabatic cooling resulting from wind expansion is 
presently neglected in our models (work under way). 
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the rate equations, assuring that the constraints of electron ther- 
ma l balance and radiat ive equilib rium are ph ysically equivalent 
rcf. lHiUier & Miller 1998 and Hillier "2003' where further dis- 
cussion concerning both methods and their correspondance is 
given). Insofar, we can use the achieved flux-conservation as 
an almost independent tool to check whether our models have 
been constructed in a consistent way. In most of the cases con- 
sidered so far we have found a perfect conservation, but in the 
worst cases (below 5% of all models) a violation up to 1.5% is 
possible. 

The third point to be discussed is mainly relevant for our 
specific approach of modeling stellar atmospheres. Presently, 
and in accordance with the majority of similar codes, we do 
not update the photospheric density stratification once it has 
been calculated. Since the photospheric structure equations are 
solved for the gas-pressure P and the density is calculated from 
the ratio P/T, the density is only as good as the initial "guess" 
for the temperature stratification. Moreover, an implicit depen- 
dence of the final temperature distribution on this initial guess 
is created. 

Thus, it is still important to obtain a fair approximation 
for the latter quantity, which in our models is accomplished 
via the corresponding NLTE Hopf-parameters (see Paper I) 
which have to account for line-blanketing effects. Meanwhile, 
we have accumulated a large set of these parameters from our 
model-grid calculations (and, for cooler temperatures, from 
corresponding Kurucz-models). If, on the other hand, the initial 
(photospheric) temperature stratification were not appropriate, 
both occupation numbers and line profiles would be aff'ected 
from the erroneous density (although the flux would be con- 
served, see above). 

In Fig.^lwe show some of our results in comparison with 
results calculated by means of WM-Basic, a code which also 
uses the electron thermal balance. Obviously, the differences 
are tiny and visible only for the temperature bumps of super- 
giants, which are predicted to be more prominent by WM-Basic. 
Note, however, that our solution is more consistent with the re- 
sults from CMFGEN (see Fig.ll4>. which will be presented in the 
next section. 

Comparing the computation time of models with and with- 
out consistent temperature structure, we find a typical differ- 
ence of a factor of two. Interestingly, the number of iterations 
becomes only moderately larger (because of the fast conver- 
gence of the temperature w hen using the electron thermal bal- 
ance, see lKubat et al."l99^). and most of the additional time is 
spent for solving the NLTE equations for the important back- 
ground elements. 

We finish this section with an interesting finding and warn- 
ing. After having calculated a large number of models with 
our code, in certain domains of Teft we have found tempera- 
ture bumps of extreme extent. In contrast to "normal" bumps 
(arising from line-heating in the outer photosphere) which are 
of the order of 2,000 K or less for O-stars (Fig. ll2> . correspond- 
ing values at lower effective temperatures might reach 5,000 K, 
as shown for an exemplary dwarf-model at Tefj - 22,400 K in 

Fig.ini 

We like to stress the fact that this behaviour has been con- 
firmed by calculations performed by means of wM-Basic and 
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Fig. 12. FASTWiND (bold) vs. WM-Basic (grey): comparison of 
temperature stratification for some of the models described in 
Sect.lsH 

Teff: 22400 K, log g: 3.90, Mdot: 1 .e-7 Msun/yr 
70 F ' ' ' ' 3 




10^ 10° 10^ 10'^ 10*^ 10® 

tau_Ross 

Fig. 13. Extreme temperature-bump around 22,000 K: fastwind 
(bold) vs. WM-Basic (grey, dashed) and cmfgen (crosses). See 
text. 

CMFGEN, kindly provided by T. Hoffmann and F. Najarro on our 
request. This finding (identical results for an unexpected and 
somewhat strange effect) allows for two conclusions. First, the 
effect is "real", at least in terms of the applied physics (see be- 
low), and second, the results from different codes using differ- 
ent techniques are strongly converging, which is very promis- 
ing and allows for an increasing trustworthiness of the results 
themselves. 

After some investigations, it turned out that the feature 
under discussion originates from bound-bound heating by 
Cm "*(which is a major ion at these temperatures), contributed 
by few transitions connected to the ground-state (singlet), to 
the meta-stable level (lowermost triplet state) and the transition 
between ground and meta-stable level at roughly 1909 A. Note 
that the latter transition has been identified to be of significant 

at this specific temperature, bumps at other temperatures originate 
from different ions, e.g., helium. 
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5V/3V 




5la/2la 





Fig. 14. As Fig. [21 but for fastwind (bold) vs. cmfgen (grey, 
dashed). The stellar parameters are similar to the models dis- 
played in Fig. [m with reff(6V) = 35861 K, reff(3V) = 
43511 K, Teff (5Ia) = 35673 K and Teff (2Ia) ^ 44642 K. cmf- 
gen re sults from the model grid as calculated bv lLenorzer et alJ 
iI20Q4 . 



importance for the energy -b alance in the wind of P Cyg, in that 
case as a cooling agent (cf. lOrewll 19851 Fig. 3). In our case, 
however, the Cm ground-state is strongly underpopulated in the 
transonic region (because of the same effect under-popu lating 
the Hen ground-state in hot stars, cf. iGabler et alJll989 ). so 
that the bracket in Eq.|92]becomes very large and the heating- 
rate enormous, also because of the large collisional strengths 
of these transitions. If, on the other hand, the contributions by 
Cm are neglected at all, a temperature bump of only moderate 
size is created. 

The lesson we learn from this exercise is two-fold. First, 
only a couple of lines (from one ion) can lead to a consider- 
able heating in stellar atmospheres, at least theoretically. Since 
this heating takes place in the outer photosphere it will have 
a significant effect on the spectra, and we can check this pre- 
diction observationally. However, we have also to consider that 
the degree of heating (i.e., the extent of the temperature bump) 
depends strongly on the corresponding collision strengths of 
the responsible transitions (as a function of temperature), and 
before relying on our results we have to carefully check for 
possible uncertainties.''' 

8. Comparison with CIVIFGEN 

In this section, we will compare the results from our models 
with corresponding results from cmfgen, with particular em- 
phasis on the optical H/He profiles which cannot be compared 
to results from WM-Basic, due to lack of comoving frame trans- 
port and adequate line-broadenin g. For this pu rpose, we have 
used the CMFGEN-simulations by iLenorzer et al. ( 2QQ4 ), who 
have provided a grid of dwarf, giant and supergiant models 
(no clumping) in the 0-/early B-star range. The correspond- 



" Note that even some of the hydrogen collision strengths have been 
revised recently, cf. .Przvbilla & Butler. 12004) . 



Fig. 15. As Fig. 1^1 but for fastwind vs. cmfgen (grey). Effective 
temperatures as in Fig. [H] Only the EUV part is plotted, at 
larger wavelength the results are extremely similar. 



ing FASTWIND models have been calculated with identical pa- 
rameters, and the explicit elements (H/He) have been treated 
with comoving frame transport. Thus, the only "physical" dif- 
ference in both calculations concerns the photospheric density 
stratification, which is approximated by a constant scale-height 
in CMFGEN, but described consistently by fastwind (cf. Sect.|5}- 

The corresponding temperature profiles are displayed in 
Fig. ^5 for two dwarf and two supergiant models with pa- 
rameters similar to our comparison with WM-Basic. Remember 
that the temperature structure is derived from radiative equilib- 
rium in CMFGEN, whereas fastwind uses the thermal balance of 
electrons in the outer atmosphere. Overall, the differences are 
small, and the extent of the temperature bumps are comparable. 
The only disagreement is found in the outer wind, where fast- 
wind uses an artificial cut-off (Tmin -OA Teff ) in order to pre- 
vent numerical problems at lower effective temperatures. We 
have convinced ourselves that this cut-off has no further conse- 
quences for the models as described here, which neglect adia- 
batic cooling in the outer wind anyway. 

Fig- El compares the corresponding EUV-fluxes, in anal- 
ogy to Fig.|8l As akeady discussed in Sect. 15.31 the largest dif- 
ferences occur in the Hen -continua. This effect can be seen 
even clearer in Fig. El lowermost panel. Regarding the super- 
giants, the deviation is contrary to our comparison with wm- 
Basic. The WM-Basic Hen -fluxes were mostly lower than those 
from FASTWIND, whereas the CMFGEN-fluxes are larger, particu- 
larly at the edges, so that the corresponding Zanstra integrals 
become larger as well. Thus, the fastwind results for Qeeii lie 
roughly in the middle of the results from cmfgen and WM-Basic, 
at least for the supergiants. Again, we like to point out the ex- 
treme sensitivity of the model predictions in this frequential 
range and warn the reader about any uncritical use of corre- 
sponding results, e.g., with respect to nebula modeling. 

Regarding the dwarf models, both codes give more or less 
identical results for the Heii -continua for Teff < 36,000 K, 
whereas at hotter temperatures extreme differences are found 
for the two models at T^s = 41,000 K and 43,500 K, respec- 
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Fig. 16. As Fig. 121 right panel, but for cmfgen vs. fastwind (pos- 
itive values result from Zanstra-integrals being larger in cm- 
fgen). Triangles: dwarfs; asterisks: supergiants. Note that the 
X-axis extends only until Teff =50,000 K. The three objects 
denoted by arrows in the lowermost panel (Hen) correspond 
to the dwarf models '3V' and '4V' and the supergiant model 
' lla', respectively. The resulting differences for A log QhcII are 
given beneath the arrows. 



tively. In contrast to both our predictions and those from wm- 
Basic the CMFOEN-models do not show any Hen -edge at all, 
cf. Fig.ini model "3V". 

Concerning the Oii-continua (actually, for the complete 
range within 300 k < A < 400 A), the hotter models (Jeff > 
35,000 K) show a higher flux-level in cmfgen, for both the su- 
pergiants and the dwarfs. We have already commented on this 
problem in Sect. 15.31 and speculated that this behaviour is re- 
lated to missing line-opacity. (On itself plays no role at these 
temperatures). Of course, we cannot exclude a problem in our 
approximate treatment of line-blocking. Finally, and in accor- 
dance with the comparison with WM-Basic, the agreement of 
the Hi - and Hei -continua is almost perfect. 



Figs. El and [121 are now the most interesting plots in 
this section, displaying the strategic H/He lines in the optical 
(cMFGEN-profiles in magenta). Regarding the dwarfs, the agree- 
ment of almost all lines is excellent. The only differences are 
found for the line cores of Hen 4686, which are shallower in 
CMFGEN at almost all temperatures, and for the Hei singlets 
for models '4V' to'6V' with Teff lying in the range between 
41,000 K and 36,000 K, respectively. (Note that for model 
'4V' Hei 4387 agrees well whereas Hei 4922 and Hei 6678 dif- 
fer). Most prominent are the differences for models '5V' and 
'6V' (the same is true for the giant models not displayed here), 
where all singlet lines predicted by cmfgen are almost a fac- 
tor of two smaller in equivalent width than those predicted by 
fastwind. Most interestingly, however, the triplet lines agree 
perfectly throughout the grid. 

So far, the origin of this discrepancy could not be identified; 
particularly, the atomic data used (incl. broadening functions) 
are very similar, and also the ionizing continua (important for 
the singlet-formation) agree very well, as shown above. One 
might speculate that there is a connection to the flux differences 
around the Hen resonance hne at 304 A or to possible discrep- 
ancies at the Hei resonance line(s), but this has to be checked 
carefully (investigations under way). Further comments on this 
discrepancy will be given after we have discussed the results 
for the supergiants. 

The corresponding profiles are displayed in Fig.^l upper 
panel. There, the situation is somewhat different from the dwarf 
case. At first, we note that the deviations of the Hei singlets 
are not as extreme as before. Significant disagreement is found 
only for Hei 4922 and 6678 (no problem for Hei 4387) in 
model "5Ia" (36,000 K), where these singlets are weak any- 
way. For model "6Ia" the differences are moderate, much less 
than the factor of two in equivalent width encountered above. 
Noticeable differences are found for other lines though. At 
first, the hydrogen Balmer line wings predicted by cmfgen are 
much stronger, which would lead to lower gravities if an anal- 
ysis of observed spectra were performed. Second, both Ho' and 
Hen 4686 show stronger wind emission which would lead to 
lower mass-loss rates compared to fastwind. Note however that 
the wind emission in both lines is a strongly increasing function 
of mass-loss (e.g.. |Puls et alJl99^ . and an analysis of observed 
spectra would result in M -differences not exceeding the 20 to 
30% level. 

The difference in the Balmer line wings points to a prob- 
lem mentioned above, namely the assumption of a constant 
photospheric scale height in cmfgen. In order to obtain an im- 
pression in how far this approximation (as well as the some- 
what artificial transition from photosphere to wind) has an in- 
fluenc e on the resulting models and profiles, iLenorzer et alJ 
(2004) have calculated an additional set of "low-gravity" su- 
pergiants, where the gravity has been lowered by typically 0. 1 
to 0.2 dex (model series "Jg") with respect to their "standard" 
grid of supergiants. Due to this manipulation, at least part of 
the effect of photospheric radiation pressure g^^i is accounted 
for (although this quantity is not constant throughout the pho- 
tosphere), since the profiles provide a measure of the effective 
gravity (i.e., ggi-^v - gmd) alone. 
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Fig. 17 . FASTWiND (black) vs . cmfgen (magenta): comparison of strategic H/He lines in the optical for the dwarf-models from the 
grid bv lLenorzer et al ]|2004). For both models, the lines have been degraded to a resolution of 10,000 and rotationally broadened 
with vsini=80 km s~'. Hei/l/l 4387, 4922 and 6678 are singlet lines, and Hei/l/l 4471 and 4713 are triplets. The horizontal and 
vertical lines in the bottom right corner indicate the scale used and correspond to 20 A in wavelength and 0.5 in units of the 
continuum, respectively (extending from 0.65 to 1.15). 
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Fig. 18. Wind-strengt h parameter Q as an optical depth invariant: H/He profiles for the model of a Cam as determined by 
[Repolust et al.' (2004), with M =6.04- lO^^Mo/yr and =32.5M0 . Overplotted in magenta are the corresponding profiles for a 
model with identical g-parameter (Eq.|93} but different mass-loss rate and radius (M =3.3-10"^MQ/yr and R^, -2\.1Mq ). 



In Fig. ^] lower panel, we compare the fastwind pro- 
files (identical to those from the upper panel, since our "high 
gravity" models do include the photospheric ^lad) with these 
low-gravity models calculated by cmfgen. Consequently, the 
photospheric densities should be much more similar than in 
the previous case, at least in those regions where the Balmer 
line wings are formed. Indeed, the differences in Hy and H^ 



have now vanished, and also the Ha emission is very simi- 
lar, except for the hottest models on the blue side of the pro- 
file. In some cases, the discrepancy for Hen 4686 has become 
weaker as well. The He: triplets have not changed (they seem 
to be almost independent on the photospheric density in cmf- 
gen), whereas a strong influence on the He: singlets is found. 
In the "critical" temperature region, they have become signifi- 




Fig. 19. As Fig. ^] however for the supergiants from the grid by iLenorzer et alJ ( l2004h . Upper panel: cmfgen models with 
"standard" gravities. Note that the differences in the wings of the Balmer lines and in the Ha -emission almost vanish if our 
results are compared to the "low gravity" cmfgen models in the lower panel (see text). 
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cantly weaker, and a strong discrepancy also for the low-gravity 
model '6IaJg' is present again, by the same degree as we have 
found for the dwarfs. 

In summary, we find a very good agreement with the op- 
tical spectra from cmfgen if the problem of different density 
stratifications is accounted for The only disturbing fact is the 
strong difference in the Hei singlets for dwarfs between 36,000 
to 41,000 K and for supergiants between 31000 to 35,000 K. 

Although it is presently not clear which profiles are "cor- 
rect" or whether the truth lies in the "middle", we like to 

S ioint out the followi ng. In our analyses of Galactic O-stars 
Repolust et alJl2004l) . no problems were found to match both 
the observed singlet and triplet lines in dwarfs. Concerning 
the supergiants, we actually met a problem for almost all stars 
cooler than 06, namely the well-known "generaliz ed dilution 
effect " (see the discussion and references in .Repolust et al.l 
I2OO4I) . Briefly, we could fit all Hei lines (singlets and triplets) 
in parallel with the Hen lines, except for Hei 4471 (triplet) 
which was predicted to be too weak. One might argue that 
this is a symptom of generally incorrect Hei lines, and specu- 
late that this problem is related to the inconsistency seen here. 
Assuming that the Hei -singlets produced by fastwind are er- 
roneous it might then be possible to fit all Hei singlets and 
the A 4471 triplet at cooler temperatures. In this case, however, 
we (and cmfgen!) would meet the problem that the other triplet 
Unes would be too strong and the Hen lines too weak! 

Presently, there is no other way out of this dilemma than to 
perform a number of detailed comparisons, with respect to both 
the models and the observations. Since the actual problem con- 
cerns the ratio of triplet to singlet lines and the problem is most 
pronounced for dwarfs, it should be possible to find a solution 
by comparing the theoretical predictions for this ratio (in terms 
of equivalent widths) as a function of T^ff vs. the observed ratio 
as a function of spectral type for a large sample of stars. Such 
work is in progress now. 

9. Model grids 

As already outlined in Sect. [2 the parameter space to be in- 
vestigated for the analysis of one object alone is large and al- 
most prohibitive for the detailed analysis of very large samples 
of stars which have recently been collected (e.g., by means of 
the multi-object spectrograph flames). Alternatively, a some- 
what coarser analysis by means of the "traditional" model-grid 
method is still applicable if an appropriate grid can be con- 
structed. In this section, we will give some suggestions for this 
objective and report on first progress. 

Although the presence of a wind introduces a large number 
of additional parameters to be considered in a fine fit (M , Voo , y6 
and R-t. ), there is a fortunate circumstance which allows for the 
construction of such model-grids with only one more parameter 
compared to grids from hydrostatic, plane-parallel models, at 
least if we do not aim at the analysis of specific (UV) resonance 
lines. 

As has been shown by, e.g.. Puis et al. (1996, see also 
ISchmutz et aPl 19891 Ide Koter et alJI 19981 for diversifications), 
the wind-emission from recombination dominated transitions 
(so-called p^-hnes) remains rather unaff'ected from the specific 



choice of the individual values of M , Vco and 7?* as long as the 
wind-strength parameter Q (also denoted as the "optical depth 
invariant"). 



Q 



M 



(Voo/?*)2 



(93) 



does not vary. In this case then, most of the other lines also 
preserve their shape. An example is given in Fig. where 
we have varied the mass-loss rate of a model of a Cam (cf. 
iRepolust et alJ i2004) by a factor of two (and accordingly the 
radius by 2'-^) without almost any effect on the resulting H/He 
spectrum. 

This behaviour (i.e., spectrum (and emergent fluxes!) de- 
pend almost exclusively on Q and not on its individual con- 
stituents) follows from the fact that 

- p^-dependent line processes (e.g., recombination lines and 
resonance lines from ions one stage below the major one^**) 
scale with Q in the wind regime, 

- the wind density scales with M/(vooR^) (continuity equa- 
tion), and 

- (res onance) lines fr om major ions scale with M/(v^7?*), 
e.g. lHam'anr] ( ll98lh . 

Thus, the common power of "1.5" with respect to Vk, and 
Ri, used in Q, which refers to the scaling of p^ lines, is also 
the best compromise to deal with the other physical parameters 
affecting a stellar model (most importantly, the Une-blocking 
which depends both on density and line opacity). 

Exploiting this knowledge, we have constructed a set of 
nine model-grids for the analysis of H/He profiles with three 
different helium abundances, Ym = 0. 1, 0. 15 and 0.2, and three 
different background metallicities, z = 1.0,0.5 and 0.2 (cf. 
Sect. 13, respectively. Each grid with given helium abundance 
and metallicity is three-dimensional with respect to the param- 
eters Tgff , log g and log Q, and the grid-spacing is roughly 
equidistant. The individual values for parameters incorporated 
into log Q (which are actually needed to calculate a specific 
model) and additional ones have been assumed according to 
present knowledge: 

- Ri, from "empirical" values, as a function of spectral type 
(Teff ) and luminosity class (logg). 

- v'co as a function of photospheric escape velocity Vesc, 



(94) 



in accordance with the results collected by 
iKudritzki&Puls'teOOOl). 

- Vel ocity exponent 13 from empirical values (see also 
Eudritzki & Pulsll200d and references therein), with /3 - 
0.9 (as a compromise) for O-stars and increasing values to- 
wards later types. 

- Mfrom log Q, and Voo as specified above. 

- micro-turbulence Vmrb = 15 km s"' throughout the grid as a 
compromise between O and B stars. 
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HEI4471: series A, log Q = -14.0000 HEI447-': series E, log Q = -12.4500 




2.0x10'^ 2.5x10* 3.0x10'^ 3. 5x10** 4.0x1o'^ 4.5x10'* S.Oxlo'^ 2.0x10'* 2.5x10* 3.0x10* 3.5x10* 4.0x10* 4.5x10* 5.0x10* 

Teff Teff 

Fig. 20. Iso-contours of equivalent widths for Hei4471, as predicted by fastwind, using results from our model-grid with Helium 
abundance Fee =0.1 and solar metallicity for background elements. Left: negligible wind. Right: typical O-supergiant wind. 
Note the effect of wind-emission, shifting the iso-contours to the left. The locations of the corresponding models are indicated 
by asterisks. 



Our present grids comprise the range 20,000 K < Tgff < 
50,000 K with A Teff = 2,500 K, loggbetween 2.2 < logg< 
3.2 at Teff = 20,000 K and 4.0 < log g < 4.5 at T^f{ = 50,000 K. 

The position of all models can be inferred from Fig. |20l 
With respect to log Q we have used values with -14.0 < log Q 
< -1 1.4 (A log Q = 0.35 in most cases), where the lowest value 
corresponds to an almost negligible wind and the highest one 
to almost Wolf-Rayet conditions. 

The denotation is such that we specify a letter for the wind 
density ("A" to "H", with densities log Q = -14.0, -13.5,-13.15,- 
12.8,-12.45,-12.1,-11.75,-11.4, respectively, if Mis calculated 
in Mo/yr, Vooin km s^'and /?* in Rq). Effective temperature 
and gravity are denoted by two ciphers each. Thus, model 
"E2730" refers to log Q = -12.45, = 27,500 K and logg = 
3.0. Typical O-type supergiants correspond to series "E", and 
typical B-type supergiants to series "D". 

For all these models we have calculated H/He profiles and 
equivalent widths in the optical and the IR. Thus, by simply 
over-plotting observed vs. simulated spectra one finds an im- 
mediate guess for the parameters Teff , log g , Fee and wind- 
strength if the background metallicity is specified and the the- 
oretical profiles have been convolved accounting for rotational 
broadening and resolution. In this way, the coarse analysis of 
one star is possible within a couple of minutes and might be 
fine-tuned by calculating specific models (particularly with re- 
spect to /3 if inferable from the emission line shapes). 

In addition, a plot of various iso-contours of calculated 
equivalent widths gives deeper insight into certain dependen- 
cies. As an example. Fig. 1201 shows the eff'ect of wind emis- 
sion on Hei4471. Further examples, particularly with respect 
to the spectral type classification criteriu m of O-star, log W = 
log(ew4471) - log(ew4541), are given in lMassev et al.l(l2005l) . 

We intend to make these grids publicly available in the near 
future when the problem regarding the Hei singlets has been 
solved. 



10. Summary and outlook 

In this paper we have described all updates applied to our previ- 
ous version of fastwind (Paper I), regarding the approximative 
treatment of metal line-blocking/blanketing and the calculation 
of a consistent temperature structure. 

The problem of line-blocking has been tackled in two steps. 
First, the occupation numbers of back-ground elements are 
calculated by an approximative solution of the corresponding 
equations of statistical equilibrium with the option that the 
most abundant elements are treated almost "exactly", i.e., by 
means of the Sobolev transport for line processes. Compared to 
alternative approaches (cf. Sect. |4} our method allows for the 
treatment of different spin systems, radially and frequency de- 
pendent radiation temperatures and a consistent ALI-iteration 
scheme. We have tested our solutions by comparing the approx- 
imative results with results from exact solutions and have not 
found any major discrepancies. 

The occupation numbers derived in this way are subse- 
quently used to calculate the Une-blocked radiation field, again, 
in an approximative way. To this end, we have formulated suit- 
able means for the opacities (in analogy to Rosseland means 
but for frequency intervals not larger than 1,000... 1,500 
km s"') and emissivities (two-level-atom approach), and the 
resulting pseudo-continuum of overlapping lines is treated by 
means of a conventional continuum radiative transfer Specific 
problems inherent in our approach (regarding a rigorous statis- 
tical description) have been pointed out and might lead to in- 
accurate solutions in a few cases. Investigations to improve our 
approach are presently under way in our group, as discussed in 
Sect. 13 

Our new version of fastwind allows for the calculation of a 
consistent temperature structure by applying a flux-correction 
method in the lower atmosphere and the electron thermal bal- 
ance in the outer one. Regarding optical H/He lines, no major 
diff'erences have been found compared to our previous NLTE 
Hopf-function method (cf. Paper I and .Repolust et al..2004i) . 
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Due to the approximations applied and as intended, the per- 
formance of our code is very fast. The total computational time 
(starting all models from scratch) is of the order of 30 minutes 
on a PC with a 2 GHz processor if only H and He lines are con- 
sider ed as exphcit ion s, whereas the inclusion of other elements 
(e.g.. lUrbaneial2004 into the "explicit" treatment requires ad- 
ditional 5 to 10 minutes each. 

The new methods have been extensively tested by compar- 
ing with results from WM-Basic and cmfgen, concerning tem- 
perature stratification, fluxes, number of ionizing photons and 
optical^' H/He profiles (comparison with cmfgen only). 

We have highlighted the importance of photospheric line- 
pressure, which is incorporated into the fastwind models and 
neglected in the standard version of cmfgen, if not coupled to 
the plane-parallel code tlusty (see Sect. Particularly, we 
have found indications that the use of the Sobolev approxima- 
tion (within the force-multiplier concept) in WM-Basic can lead 
to an underestimate o f this quantity, as already predicted by 
IPauldrach et alJ(ll986l) . On the other hand, the density/velocity 
stratification resulting from our approach (smoothly connecting 
the quasi-static photosphere and a y6-law wind) agrees surpris- 
ingly well with the hydrodynamic structure as calculated from 
a consistent solution if yS is not too different from the "canoni- 
cal" value of 0.8. . . 1.0. 

All three codes predict almost identical temperature struc- 
tures and fluxes for A > 400 A, whereas at lower wavelengths 
certain discrepancies are found. Compared to WM-Basic (us- 
ing an identical line list for the background elements), our 
supergiant models differ only in the Heii continua, where 
the FASTWiND-fluxes are somewhat larger, but still lower than 
the corresponding fluxes from cmfgen. Since fluxes and corre- 
sponding numbers of ionizing photons react extremely sensi- 
tive to subtle model differences in this wavelength regime, we 
consider any uncritical use of these quantities as being danger- 
ous. 

Major discrepancies are also found in the range 300 A < A 
< 400 A, i.e., in the On continuum and at the Hen 304 reso- 
nance line. Compared to both WM-Basic and cmfgen, our dwarf 
models produce less flux in this region (more blocking or less 
re-emission), whereas the supergiant models of fastwind and 
WM-Basic agree very well. The supergiant models of cmfgen, 
on the other hand, show much less blocking which might point 
to some missing opacity. Again, we like to stress that the Hi 
and Hei continua agree very well in all three codes. 

For the optical H/He lines, the coincidence between fast- 
wind and CMFGEN is remarkable, except for the Hei singlets in 
the temperature range between 36,000 to 41,000 K for dwarfs 
and between 31,000 to 35,000 K for supergiants, where cmfgen 
predicts much weaker singlets. Up to now, the origin of this 
discrepancy could not be identified, but work is under way to 
solve this problem. 

Although it is reassuring that the different codes agree 
well with respect to most of their predictions, this is only 
part of the story. One particularly disturbing fact concerns the 



^' IR-lines will be presented in a forthcoming paper (Repolust et al., 
in prep, for A&A), with a similar agreement between fastwind and 
CMFGEN as for the optical ones. 



present mismatch between the parameters obtained from an 
analysis in the optical and the UV, respectively. In the major- 
ity of cases, the UV gives lower effective temperatures, i.e., 
of the order of 2,000 to 4,000 K, if one compares th e anal- 
yse s of Galactic stars perfor med by l^ianchi & Garcia | (|2002| ) 
and ' Garcia & Bianchil ( l2004l) with results from Repolust et alJ 
(I2004ir (wM-Basic vs. fastwind) an d the corresp o nding work 
for MageUanic C loud stars by iHillier et all (l2003l) and 
iBouret et al.1 ll2003ll ( partly including also the optical range) 
with the results from iMassev et all 1 2QQ4i |2 QQ^) (cmf gen vs. 
fastwind). (Interestingly, the work by Crowther et alJ20 02 (cm- 
fgen) indi cates high er temperatures for MC supergiants than 
derived bv lMassevet al„2005. ) 

Part of this discrepancy (if combined UV/optical analyses 
are compared) might be related to the Hei singlet vs. triplet 
problem as discussed above. Note, however, that this would ac- 
count only for discrepancies in certain domains of the Tetf space 
and would typically result in maximum differences of the or- 
der of 2,000 K, as has been found from a number of test cal- 
culations performed by one of us (J.P) and F. Najarro (using 
cmfgen), which will be reported on in a forthcoming publica- 
tion. Mo reover, the temperatu re scale for O-type dwarfs as de- 
rived by iMartins et alJ ( l2002l) using cmfgen and concentrating 
on the classification criterium Hei 4471 (triplet) vs. Heii4541 
is actually 1,000 to 2,000 K hotter than the calibration by 
iRenohistet al. (20o4 . 

In a recent paper. iMartins et alJ (l2004t have discussed the 
uncertainties in Teff which is obtained by relying on different 
diagnostic tools in the UV, analyzing four SMC-N81 dwarfs 
of spectral types 06.5 to 08.5. From the specific values de- 
rived from the UV-color index, the ionization balance of Oiv/v 
and Feiv/y and flie Nvl238/1242 and Cral426/1428 doublets, 
respectively, they quote a typical uncertainty of + 3,000 K in 
Teff , which might easily account for part of the discrepancies 
with the optical. 

Unfortunately, it is rather difficult to compare the differ- 
ences obtained so far in a strict one-to-one case, simply be- 
cause the corresponding samples hardly overlap. In particu- 
lar, a large fraction of the objects which have been analyzed 
by means of cmfgen are somewhat ex treme, comprising eit her 
supergiants with (very) dens e winds jCrowther et alJl2002h or 
dwarfs with very thin winds (Martins et al.''2004V The analysis 
of SMC stars by .Bouret et al. 12003 ). on the other hand, cov- 
ers only a sample of 6 dwarfs, in contrast to the larger sample 

I ' "-I I 1 1 — - — I 

by Massey et al. (2004, 2005), and, therefore, it is not clear in 

how far selection effects do play a role. Finally, it is interesting 
to note that at least for one object in common, the 04I(f) star ^ 
Pup (HD 668 11), the different analyses give almost identical re- 
sults ( Crowther et al 2002, Repolust et al. 2 004 and Pauldrach 

et al., in prep, for A&A, analyzing the UV by means of wm- 
Basic). 

Thus, we conclude that the present status of hot star pa- 
rameters is not as clear as we would like it to be. Actually, we 
need to understand a number of additional physical processes 
and their influence on the derived parameters. Most important 
are the direct and indirect effects of the Une-driven wind insta- 
bility, i.e., the formation and interaction of clumps and shocks 
leading to X-ray emission and enhanced EUV-flux in the wind 
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(e.g., lFeldmeieret"aDll997t IPauldrach etaPbOOll) . Although 
incorporated to some extent into present codes, there are sim- 
ply too many questions to be answered before we can consider 
these problems as solved. To give only two examples: We do 
not know the spatial distribution of the "clumping factor", and 
also the X- ray emission is only on the verge of bei ng under- 
stood (e.g..' Kramer et alJ2003UOskinova et alJ2004 . 

Before these effects can be treated in a realistic way, we 
suggest to primarily rely on diagnostic tools which are least 
"contaminated", i.e., to concentrate on weak lines formed in the 
stellar photospheres (except, of course, the mass-loss indicators 
which will always be affected by clumping). Future investiga- 
tions of O-type stars performed by fastwind will have to utilize 
not only H and He but also metal lines, as already incorporated 
into the analysis of B-stars (cf. SectQ. Particularly, one of the 
most important tools will be nitrogen with its strong sensitivity 
even at higher temperatures where Hei begins to fail. Work in 
this direction is under way. 
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